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SECTION  I 


INTRODUCTION 


The  major  goal  o£  the  research  described  In  this  report  is  the  calculation 
of  electron- impact  dissociation  cross  sections  for  vibrationally  and  rotation- 
ally  excited  diatomic  molecules.  Dissociation  by  electron  impact  to  form 
neutral  species  is  an  important  process  in  gas-phase  electric  discharges  and 
can  significantly  affect  the  efficiency  of  laser  plasmas.  Although  large 
populations  of  rotationally  and  vibrationally  excited  species  exist  in  plasmas, 
little  is  known  either  experimentally  or  theoretically  about  the  dependence 
of  electron-impact  dissociation  cross  sections  upon  the  initial  rovibrational 
state  of  the  diatomic  molecule.  The  cross  sections  reported  here  are  important 
to  the  basic  understanding  of  the  dissociation  process  as  well  as  for  input 
to  modelling  studies  of  gas-phase  electric  discharges. 

The  first  task  was  the  selection  and  verification  of  a  theoretical  method 
that  was  capable  of  reasonably  accurate  predictions  of  cross  sections  for 
dissociation  by  electron  impact  at  relatively  low  collision  energies.  Several 
methods  have  been  previously  applied  to  electron-impact  excitation  and  dissocia¬ 
tion,  but  none  of  these  have  studied  the  effects  of  changing  the  initial  rota¬ 
tional  and  vibrational  state  upon  the  cross  sections.  From  these  methods, 
the  impact  parameter  method  was  selected  as  the  one  best  suited  to  our  needs. 
We  extended  it  to  treat  dissociation  and  to  include  the  effects  of  vibrational 
and  rotational  motion.  The  method  was  verified  by  application  to  transitions 
in  the  H2  molecule  for  which  experimental  results  were  available  for  comparison. 

The  impact  parameter  method  was  then  applied  to  transitions  in  the  O2, 
S2,  HC1,  and  HgBr  molecules  for  electron  translational  energies  in  the  range 
from  threshold  to  25  eV.  First,  a  survey  was  made  of  the  dependence  of  the 
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cross  sections  on  initial  vibrational  state.  For  these  studies  the  rotational 
degrees  of  freedom  were  treated  as  being  degenerate.  Cross  sections  were 
computed  for  Initial  vibrational  states  “  0*  1*  and  2  for  O^and  HC1  and 
for  m  0  to  4  for  and  HgBr.  Then  a  survey  was  made  of  the  dependence  of 
the  cross  sections  on  initial  rotational  state  for  various  initial  vibrational 
states.  Cross  sections  were  computed  for  a  distribution  of  initial  rotational 
states  characterized  by  temperatures  of  0,  300,  400,  600,  and  1000  K. 

The  results  of  the  calculations  are  summarized  in  Section  II.  This 
section  also  provides  a  summary  of  the  general  trends  observed  in  our  studies 
and  assesses  the  general  applicability  of  the  calculations.  The  Appendices 
contains  manuscripts  prepared  for  publication,  which  give  more  detailed 
descriptions  of  the  methods  employed,  the  calculations,  and  the  results. 
Appendix  A  describes  the  impact  parameter  method,  our  extensions  to  it,  and 
the  application  of  the  method  to  electronic  transitions  in  the  H2  molecule. 
Appendix  B  describes  the  electronic  structure  calculations  for  0^  and  S2»  and 
Appendix  C  presents  the  dynamical  calculations  on  these  molecules.  Appendix  D 
provides  a  test  of  some  of  the  methods  by  comparing  our  calculation  of 
oscillator  strengths  with  experimental  ones  for  the  well-studied 
Schumann-Runge  transition  in  0^.  Appendix  E  and  F  present  our  calculations  on 
the  HC1  and  HgBr  molecules,  respectively.  Finally,  flowcharts  of  the  computer 
programs  developed  for  the  impact  parameter  method  are  provided. 


SECTION  II 


SUMMARY  OF  WORK  PERFORMED 

1.  The  Impact  Parameter  Method 

The  impact  parameter  (IP)  method  for  diatomic  molecules  as  originally 
formulated  by  Nazi  is  extended  to  treat  dissociation  and  to  include  the  effectB 
of  vibrational  and  rotational  excitation.  In  its  extended  form,  the  expression 
for  the  state-to-state  cross  section  retains  the  simple  form  obtained  by  Hazi; 
the  cross  section  is  the  product  of  a  structural  factor  S^f,  which  depends 
upon  the  electronic  properties  of  the  diatomic  molecules,  and  a  dynamical 
factor  D^fCE),  which  depends  upon  the  translational  energy  E  of  the  incident 
electron  and  the  transition  energy  E.  This  expression  has  a  form  similar 
to  that  for  the  photodissociation  cross  section  except  that  the  dynamical 
factor  replaces  the  factor  depending  on  the  frequency  of  the  radiation  and 
the  symmetry  dependence  in  the  structural  factor  is  different.  In  both  cases 
the  structural  factor  is  proportional  to  the  absolute  oscillator  strength, 
which  can  be  obtained  directly  from  spectroscopic  data  or  from  calculations 
based  upon  potential  energy  curves  and  the  transition  dipole  moment.  The 
dynamical  factor  is  easily  evaluated  in  terms  of  modified  Struve  and  Bessel 
functions . 

In  the  IP  method,  nuclear  motion  is  treated  within  the  adiabatic 
approximation  and  the  motion  of  the  incident  electron  is  treated  as  a 
straight-line  classical  trajectory.  The  probability  of  an  electronic  transition 
for  a  fixed  nuclear  geometry  and  initial  impact  parameter  is  given  by 
time-dependent  perturbation  theory,  where  the  molecular  electrons  are  perturbed 
by  the  incident  electron.  In  this  method  only  the  repulsive  interactions 


treated;  exchange  interactions  are  neglected.  The  coulomb  repulsion  term 
expanded  in  a  multipole  series  and  approximated  by  its  asymptotic  form. 

»  assumes  that  the  major  contributions  to  the  cross  section  come  from  large 
ict  parameters.  Cross  sections  are  obtained  by  averaging  the  probabilities 
r  molecular  orientations ,  nuc)aar  separation,  and  initial  impact  parameters, 
onzero  value  bp  is  used  as  the  lower  limit  in  the  integration  over  impact 
ameter  to  avoid  divergence  of  the  cross  section.  The  minimum  impact 
ameter  is  fitted  so  that  the  IP  cross  sections  match  those  of  the  Born 
roximation  at  very  high  translational  energies. 

Thus  the  IP  method  is  a  semiempirical  theory  designed  to  have  the  correct 
h  energy  behavior.  Hazi  and  coworkers  have  shown  that  it  also  gives  a 
ter  prediction  of  cross  sections  at  low  energies  than  plane-wave  methods 
h  as  the  Born  approximation.  Because  of  the  neglect  of  electronic  exchange 
ween  the  incident  and  molecular  electrons,  the  IP  method  is  limited  to 
dies  of  electronic  transitions  which  are  optically  allowed  and  it  becomes 
s  accurate  at  threshold  energies. 

In  principle,  the  IP  cross  sections  can  be  evaluated  using  experimental 
a  to  obtain  the  structural  factor.  However,  in  most  of  the  cases  that 
studied,  this  data  is  incomplete  and  does  not  include  information  about 
ational-state  dependence.  Therefore,  all  of  our  structural  factors  are 
culated  from  potential  energy  curves  and  dipole  transition  moment  curves. 

best  available  potential  data  are  used:  RKR  (Rydberg-Klein-Rees)  points 
n  available,  otherwise  the  most  accurate  ab  initio  points.  The  transition 
ole  moments  are  taken  from  ab  initio  calculations  since  they  are  hard  to 
ract  from  experimental  data.  As  a  check  on  the  accuracy  of  this  input 
a,  we  compare  our  computed  oscillator  strengths  with  experimental  ones 
n  they  are  available. 
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le  I.  Test  of  the  Treatment  of  Small  Impact  Parameters  in  the  IP  Method. 
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the  bg  is  adjusted  to  make  o  (E)  agree  with  the  BA  cross  section  at  high 


energies.  A  measure  of  the  importance  of  the  region  between  b  =  0  and  bg 


“IP  ~xp  ~ ip  IP 

is  then  given  by  the  ratio  of  o,  to  o  and  the  ratio  between  o,  and  a 

b0  b0 


iP, 


gives  an  estimate  of  the  error  in  o  (E)  from  the  method  of  treating  low  impact 
parameters.  In  Table  I  we  present  these  two  ratios  at  three  energies  for 
several  of  the  electronic  transitions  considered  in  this  report.  Also  included 
in  the  table  are  the  two  minimum  impact  parameters  bg  and  bg  and  the  values 


of  the  opacity  functions  at  these  values  of  the  impact  parameter.  The  largest 


IP  ~  IP 

difference  between  a  and  o.  is  for  the  X  to  II  transition  in  0o»  but  they 

b0 


IP 


differ  by  only  80%.  The  contribution  to  o  from  low  impact  parameters  is 
typically  between  207,  and  60%,  the  lower  values  associated  with  small  bg'8  • 
The  first-order  perturbation  approximation  is  most  valid  when  the  opacity 
functions  are  much  less  than  one.  Thus  the  IP  method  should  be  better 


justified  for  systems  in  which  P(bg,E)  is  small  such  as  the  X  to  IT  transition 


.3  - 


in  O2,  the  X  to  2  transition  in  S2>  the  X  to  A  transition  in  HCl,  and  the 


X  to  A  transition  in  HgBr.  The  X  to  II  transition  in  S2  is  clearly  a  case 
where  the  approximations  break  down. 

The  energetic  threshold  energy  for  dissociation  (the  lowest  energy  at 
which  dissociation  is  allowed)  can  be  accurately  computed  to  within  0.01  to 
0.1  eV  for  all  but  the  HgBr  system  using  known  spectroscopic  data.  The 
dynamical  threshold  energy  is  defined  in  this  work  as  the  lowest  energy  at 
which  the  cross  section  becomes  0.01  times  its  maximum  value.  The  dynamical 
threshold  is  determined  mostly  by  the  lowest  transition  energy  at  which  the 
overlap  between  the  initial  and  final  nuclear  wave funct ions  becomes  appreciable. 


I 
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unction  is  interpreted  as  the  probability  of  electronic  excitation  for  a 
ranslational  energy  E  and  impact  parameter  b.  The  integral  cross  section 
s  related  to  the  opacity  function  by 

o(E)  “  2fl  db  b  P(b,E)  (1) 

0 

nd  the  opacity  function  is  bounded  by  one.  However,  in  the  IP  method  P(b,E) 
liverges  for  small  b.  To  avoid  divergence  of  the  integral  cross  section, 
n  the  p  esent  work  the  lower  limit  of  integration  in  equation  (1)  is  replaced 

»y  b0 


oIP(E)  =*  2fl  J°°  db  b  PIP(b,E)  (2) 

b0 

Contributions  to  o(E)  are  excluded  from  regions  where  the  approximations  of 

the  theory  break  down  and  where  exchange  effects  become  important.  bp  is 

thosen  so  that  the  cross  section  from  the  IP  method  is  the  same  as  that  from 

the  BA  method  at  very  high  energy.  To  test  the  effects  of  this  method  for 

IP 

•electing  bp  and  to  estimate  the  possible  contribution  to  o  (E)  from  b  = 
)  to  bg  an  alternative  treatment  of  the  low  b  region  is  tested.  Instead  of 
ising  equation  (2)  which  assumes  the  contributions  to  (E)  from  b  between  0 
ind  bg  are  negligible,  the  probability  is  assumed  to  be  constant  for  this 
range  of  impact  parameters  and  the  cross  section  is  approximated  by 

oIP(E)  =  o*P  +  2n  J"  db  b  PIP(b,E)  (3) 

°  b0 

rhere 
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5. 


Error  Estimates 


It  is  difficult  to  obtain  quantitative  estimates  of  the  probable  errors 
in  the  IP  method.  The  method  contains  several  approximations  and  the  validity 
of  each  of  these  approximations  is  dependent  upon  the  specific  electronic 
transition  being  studied.  To  help  understand  the  order  of  magnitude  of  the 
errors  expected  we  compare  the  IP  method  with  the  Born  approximation  for  which 
the  limitations  are  better  understood. 

Both  methods  separate  nuclear  and  electronic  motion  by  the  Born-Oppenheimer 

t 

(fixed-nuclei)  approximation  and  this  approximation  is  very  well  justified. 

Both  are  first-order  perturbation  treatments  and,  therefore,  assume  that 

probabilities  for  electronic  excitation  for  fixed  impact  parameter  (or  orbital 

angular  momentum)  are  much  less  than  unity.  Furthermore,  both  methods  only 

treat  the  perturbation  potential  as  a  coulombic  repulsion  term  and  neglect 

the  electron  exchange  interaction.  As  pointed  out  by  Inokuti^  these 

approximations  are  valid  when  the  translational  energy  of  the  incident  electron 

is  much  larger  than  the  orbital  energies  of  the  molecular  electrons. 

Approximating  the  orbital  energies  by  the  lowest  ionization  potential  for 

the  diatomic  molecules  studied  here,  the  translational  energy  must  be  much 

greater  than  from  9  to  13  eV.  Although  this  condition  is  not  met  for  energies 

below  100-1000  eV,  the  BA  and  IP  methods  can  predict  accurate  integral  cross 

sections  to  much  lower  energies.  Typically  BA  cross  sections  are  factors 

2 

of  2  to  5  times  too  large  in  the  medium  energy  range  below  50  eV . 

The  limit  at  which  the  above  validity  criterion  is  met  can  also  be  reached 
by  increasing  the  orbital  angular  momentum  (i)  or  the  impact  parameter  (b). 
Thus  at  large  l  or  b,  the  BA  and  IP  expressions  for  the  opacity  function  P(b,E) 
and  for  the  large-angle  differential  cross  section  become  valid.  The  opacity 


section  is  seen  as  the  rotational  temperature  is  raised  from  0  to  finite 
temperatures.  This  is  a  symmetry  effect  resulting  from  the  lack  of  a  j  ■ 
0  state  for  electronic  states. 

The  impact  parameter  method  as  described  in  this  report  is  applicable 
to  any  optically  allowed  transition  in  diatomic  molecules.  The  applications 
presented  here  have  been  to  dipole  allowed  transitions,  but  the  methods  can 
be  extended  to  treat  higher  multipole  transitions. 


4. 


General  Trends 


The  cross  sections  are  most  strongly  Influenced  by  the  overlap  of  the 
Initial  vibrational  wavefunction  with  the  final  wavefunction  for  nuclear  motion 
and  the  magnitude  of  the  transition  dipole  moment.  Transitions  with  poor 
Franck-Condon  overlaps  for  ”  0  will  generally  have  a  large  enhancement 
with  increasing  vj  whereas  transitions  with  large  Franck-Condon  overlaps  for 
vj  e  0  will  generally  have  little  dependence  upon  or  will  decrease  with 
increasing  vj. 

In  general,  the  effect  of  changing  initial  rotational  states  is  much 
smaller  than  the  effect  of  increasing  the  vibrational  state.  Within  the  IP 
method,  only  transitions  in  which  the  rotational  quantum  number  j  changes 
by  ±1  are  allowed.  Therefore,  increasing  does  not  change  the  transition 
energy  appreciably  and  the  threshold  energy  for  excitation  changes  only 
slightly.  The  main  effect  of  increasing  is  to  increase  the  effective  poten¬ 
tial  by  the  centrifugal  term  J(J+l)/2pR^.  The  increase  in  the  potential  ia 
large  for  smaller  R,  so  the  potential  well  and  the  vibrational  wavefunction 
for  a  given  level  V£  are  shifted  to  large  internuclear  distances  and  changed 
in  shape  slightly.  However,  changing  is  a  small  perturbation  on  the  wave- 
function  compared  to  changing  v^. 

The  two  exceptions  to  this  behavior  that  we  observed  were  for  the  X-E 
transition  in  O2,  the  X-23j)“  transition  in  S2»  and  X-n  transitions  in  general. 
For  small  the  X-E  excitation  of  O2  is  predominately  to  quasibound  states 
of  the  local  well;  however,  as  (and  thereby  jf)  are  increased  the  local 
well  in  the  E-state  potential  becomes  more  shallow.  Eventually  it  disappears 
whereupon  the  excitation  becomes  dissociative.  Thus  the  dissociative  cross 
sections  are  enhanced  by  a  factor  of  2  as  the  rotational  temperature  is  raised 
from  0  to  300  K.  For  X-n  transitions  a  factor  of  2  decrease  in  the  cross 


The  X-to-A  dissociation  cross  sections  have  an  interesting  dependence 
upon  rotational  temperature  Trot.  Upon  increasing  Trot  from  0  to  300  K  the 
cross  sections  decrease  by  about  a  factor  of  2.  Increasing  Trot  to  1000  K 
decreases  the  cross  sections  further,  but  by  less  than  107..  This  behavior 
was  also  observed  in  calculations  for  v  ■  0,  1,  and  2. 

D.  HgBr 

2  2+  2 

The  All,  BE  ,  and  C  II  states  of  HgBr  are  the  three  lowest  excited  elec- 

2  + 

tronic  states  which  have  optically  allowed  transitions  to  the  X  E  electronic 
state.  The  main  route  for  electron- impact  dissociation  of  the  X  state  is 
through  the  A  state;  for  all  considered  the  contributions  to  the.  dissociation 
cross  section  from  the  B  and  C  states  are  less  than  10”  cm2 .  The  X-to-A 
dissociation  cross  sections  have  very  little  dependence  upon  the  initial  vibra¬ 
tional  state  V£.  Therefore,  the  overall  dissociation  cross  sections  will 
not  be  vibrationally  enhanced. 

The  X-to-A  dissociation  cross  sections  have  the  same  rotational- temperature 
dependence  as  the  X-to-A  transition  in  HCl.  As  Trot  is  increased  from  0  to 
300  K  the  cross  sections  decrease  by  about  a  factor  of  2,  and  as  Trot  is  in¬ 
creased  to  1000  K  the  cross  sections  decrease  further,  but  by  less  than  1%. 
We  also  observe  this  behavior  for  v  *  3  and  4.  Thus  we  see  a  negligible 
dependence  of  the  dissociation  cross  section  on  temperature  for  this  system. 


with  rotational  temperature,  particularly  for  v^>0.  The  X-B  dissociation 
cross  section  for  V£  =  0  decreases  by  257.  as  the  rotational  temperature  is 
increased  from  0  K  to  1000  K.  For  v^  *  1  the  cross  section  increases  by  5% 
over  this  temperature  range,  while  for  v^  ■  2  and  v^  =  3  the  increase  is  117.. 
For  vj  *  4  the  increase  is  8%.  The  X-H  cross  section  has  a  negligible 

temperature  dependence  regardless  of  initial  vibrational  state. 

3  - 

The  X-2  E  cross  section  for  v-i  *  0  increases  by  a  factor  of  1.6  from 

u  1 

0  K  to  300  K,  with  negligible  change  for  further  increase  in  temperature. 

For  *  3  the  0  K  to  300  K  change  is  a  factor  of  4,  with  the  cross  section 
then  showing  a  decrease  of  15%  by  1000  K.  For  v^  ■  4,  the  0  K  to  300  K  increase 
is  by  a  factor  of  8,  with  a  subsequent  decrease  by  26%  as  the  temperature 
is  raised  to  1000  K. 

C.  HCl 

The  main  route  for  electron-impact  dissociation  of  the  X  E  state  of 
HCl  is  through  the  A1!!  state.  The  A  state  is  the  lowest  state  accessible 

by  an  optical  transition  from  the  ground  electronic  state.  The  next  two  states 
allowed  are  the  2^E+  and  C^Jl  states.  Each  has  a  broad  well  so  that  the  energy 
at  the  ground-state  equilibrium  geometry  is  below  the  respective  dissociation 
limits.  Therefore,  excitation  to  bound  levels  will  be  preferred  over  dissocia¬ 
tion.  Furthermore,  since  these  dissociation  limits  are  both  approximately 
15  eV  above  the  X-state  minimum,  the  threshold  for  dissociative  excitation 
through  these  states  will  be  much  higher  than  that  of  the  A  state. 

The  X-to-A  dissociation  cross  sections  are  slightly  enhanced  near  threshold 
and  for  energies  above  10  eV  with  increasing  v^.  From  12  to  20  eV  the  V£ 
■  1  cross  sections  are  about  15%  larger  than  those  for  *  0,  and  the  vj 
M  2  cross  sections  are  about  10%  higher  than  those  for  *  1. 


B.  S2 


$2  is  isoelectronic  with  O2  so  the  three  analogous  excited  electronic 

3-3  3  - 

states  were  considered  for  this  system:  B  E  ,  B1  n  ,  and  2  E  .  The  avoided 

J  u  u  u 

3  - 

crossings  seen  in  the  O2  system  are  not  present  here  except  for  the  2  E^  state, 

which  has  a  local  minimum  in  its  potential  energy  curve.  Cross  sections  for 

3  - 

dissociation  of  the  X  E^  electronic  state  through  the  B  state  are  2  to  4  orders 

of  magnitude  smaller  than  those  for  transitions  to  bound  vibrational  states. 

Increasing  the  initial  vibrational  level  greatly  enhances  the  cross  sections 
for  dissociation  but  has  little  effect  upon  the  bound-state  cross  sections. 
For  example,  at  20  eV,  changing  from  0  to  1  decreases  the  bound-state  cross 
section  by  only  a  few  percent  and  increases  the  dissociative  cross  section 
by  a  factor  of  14.7.  Increasing  up  to  4  increases  the  dissociation  cross 
section  by  three  orders  of  magnitude  relative  to  =  0. 

The  II  state  is  repulsive,  and  for  vj  **  0  the  X-to-Il  transition  is  the 
dominant  route  for  dissociation.  This  cross  section  reaches  a  peak  of  4.6x10"^^ 
cm^  at  an  energy  of  6.5  eV  and  falls  off  rapidly  to  a  value  of  2x10"^^  cm^ 
at  20  eV.  The  cross  sections  are  enhanced  slightly  upon  increasing  v^;  at 
6.5  eV  the  increase  is  only  257.  in  going  from  “  0  to  4. 

3  - 

The  contributions  to  the  dissociative  cross  section  from  the  X-to-2  Eu 

transition  are  harder  to  assess  because  of  the  existance  of  predissociating 

3  - 

quasibound  states  in  the  local  well  of  the  2  E^  state.  Although  the  dissocia¬ 
tion  cross  sections  for  this  transition  show  a  large  enhancement  with  increasing 
vj  (over  4  orders  of  magnitude  from  vj  *  0  to  4),  it  will  still  be  less  than 
about  3x10”^  cm^  for  vj.  *  4. 

The  dependence  of  the  dissociation  cross  section  on  rotational  temperature 

for  S2  is  similar  to  O2  in  that  the  X-B  and  X-fl  transitions  show  little  effect, 

3  - 

regardless  of  vj,  while  the  2  E^  cross  section  shows  a  significant  increase 
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for  V£  *  0  or  1,  but  it  should  provide  between  20  to  30%  of  the  dissociative 
cross  section  for  v^  =  2  for  energies  above  about  12  eV.  For  v^  “  0  and  1 
the  total  dissociative  cross  section  rises  from  a  threshold  near  7  eV  to  about 
8x10"17  cm2  at  20  eV,  and  for  vj  ■  2  it  rises  to  maximum  of  almost  11x10“^^ 
cm^  at  20  eV. 

No  appreciable  change  is  predicted  for  the  X-to-B  dissociation  cross 
section  for  =  0,  1,  or  2  as  the  rotational  temperature  is  increased  from 

0  to  1000  K.  The  largest  changes  (about  +67.)  occur  for  near-threshold  energies, 
and  at  20  eV  the  changes  are  less  than  17..  The  X-to-B  bound-state  cross 
sections  decrease  as  the  rotational  temperature  is  increased;  typically  by 
7-8%  as  Trot  is  increased  from  0  to  300  K  and  by  25%  as  it  is  increased  from 
0  to  1000  K. 

For  the  X-n  transition  the  threshold  V£  *  0  cross  section  increases  by 
47.  from  0  K  to  300  K,  with  no  further  change  to  1000  K.  For  v^  *  1,  the  cross 
section  decreases  by  4%  between  0  K  and  300  K,  and  decreases  an  additional 
2%  by  1000  K.  At  higher  energies  the  decrease  between  0  K  and  300  K  is  7%, 
with  an  additional  2%  decrease  by  1000  K.  Thus  the  cross  section  for  this 
transition  is  not  very  sensitive  to  rotational  temperature. 

The  X-to-E  dissociation  cross  sections  show  the  largest  change  as  the 
rotational  temperature  is  increased.  At  energies  near  threshold  the  vj  * 
0  cross  section  increases  by  factors  of  2  and  2.1  upon  increasing  the 
temperature  from  0  to  300  K,  and  from  300  K  to  1000  K,  respectively.  At  higher 
energies  the  increases  are  factors  of  1.65  and  1.75  for  the  same  changes  in 
temperature.  The  v^  =  1  threshold  cross  section  decreases  by  0.8  in  going 
from  0  K  to  300  K,  with  no  further  change  observed  in  increasing  the  temperature 
to  1000  K.  At  higher  energies  the  behavior  is  similar,  except  that  the  factor 
reduces  to  0.6.  For  *  2,  the  threshold  cross  section  increases  by  1.3 
between  0  K  and  300  K,  with  no  further  change  in  going  to  higher  temperature. 
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3. 


Dynamical  Calculations 

A.  O2 

Three  excited  electronic  states  were  considered  for  this  systems  b"*E~ 

u’ 

3  - 

1  Ity,  and  E  E^.  The  presence  of  a  (diabatic)  Rydberg  state  of  the  same  symmetry 

causes  the  potential  curves  for  each  of  these  states  to  have  shoulders  or 

local  minima  resulting  from  one  or  more  avoided  crossings  near  the  equilibrium 

geometry  of  the  ground  electronic  state.  Avoided  crossings  also  cause  rapid 

changes  in  the  transition  dipole  moment  and  affect  the  shapes  of  the  excitation 

and  dissociation  cross  sections  curves  as  a  function  of  translational  energy. 

Cross  sections  for  dissociative  transitions  from  the  ground  electronic 

3  - 

state  X  E  through  the  B  state  are  2  to  4  orders  of  magnitude  larger  than 
6 

those  for  transitions  to  bound  vibrational  levels  of  the  B  state.  Increasing 
the  initial  vibrational  level  greatly  enhances  the  cross  sections  for  excitation 
to  bound  states  but  has  lLttle  effect  upon  the  dissociative  cross  sections. 
For  example,  at  20  eV,  changing  from  0  to  1  increases  the  bound-state  cross 
section  by  a  factor  of  10.1  and  increasing  vj  from  1  to  2  increases  it  by 
a  factor  of  4.2.  At  the  same  energy  the  dissociative  cross  section  decreases 
by  only  19%  upon  increasing  from  0  to  2. 

The  n  state  is  mostly  repulsive,  and  the  dissociative  cross  sections 
for  the  X-to-n  transition  are  50  to  70  times  smaller  than  the  dissociation 
cross  sections  for  X-to-B  transitions.  The  X-to-fl  cross  sections  show  only 
a  moderate  enhancement  with  increasing  v^. 

The  contributions  to  the  dissociative  cross  section  from  the  X-to-E  transi¬ 
tion  are  harder  to  assess  because  of  the  existence  of  predissociating  quasibound 
states  in  the  local  well  of  the  E  state.  Our  estimate  is  that  the  E  state 
will  not  contribute  significantly  to  the  total  dissociative  cross  section 


be  higher,  but  the  various  states  can  be  treated  equally  since  the  orbitals  can 
be  optimized  for  each  state.  Because  they  are  optimized  especially  for  the 
configuration  space  used,  the  potential  energy  curves  are  smoother  and  probably 


more  accurate. 


Ab  Initio  Electronic  Structure  Calculations 


Ab  initio  electronic  structure  calculations  are  performed  to  provide 
input  for  the  impact  parameter  calculations.  Although  RKR  potential  data 
are  available  for  many  of  the  s  ates  of  interest,  they  do  not  provide  global 
descriptions  of  the  curves.  In  particular,  the  repulsive  parts  of  the 
potentials,  which  are  critical  for  calculations  of  dissociative  cross  sections, 
must  be  obtained  from  ab  initio  calculations.  Although  the  transition  dipole 
moment  for  the  equilibrium  geometry  can  be  estimated  from  spectroscopic  data 
by  using  the  Franck-Condon  approximation,  accurate  calculations  of  oscillator 
strengths  and  electron-impact  cross  sections  require  the  transition  dipole 
moment  as  a  function  of  R.  Therefore,  it  must  be  obtained  from  ab  initio  calcu¬ 
lations.  Finally,  the  electronic  wavefunctions  and  transition  dipole  density 
matrix  are  required  as  input  to  the  Born  calculations  used  in  determining 
bg.  These  are  also  provided  by  ab  initio  calculations. 

We  have  used  the  most  accurate  ab  initio  data  available.  Some  data  has 
been  taken  from  the  literature  to  augment  or  replace  our  own  calculations. 
The  method  we  have  employed  is  the  mul-iconfigurational  self-consistent  field 
(MCSCF)  method.  It  is  well  suited  for  the  description  of  ground  and  excited 
states  of  diatomic  molecules.  MCSCF  is  a  state-of-the-art  method  "beyond" 
Hartree-Fock;  i.e.,  it  treats  electron  correlation.  The  method  not  only  opti¬ 
mizes  the  eigenvectors  (i.e.,  the  configuration  coefficients)  as  does  the 
configuration  interaction  (Cl)  method,  it  also  optimizes  the  orbitals  from 
which  the  configurations  are  constructed.  Thus  for  a  given  basis  set  and 
configuration  space,  it  provides  the  optimum  result.  Although,  for  a  given 
effort  in  computer  time,  the  configuration  space  of  an  MCSCF  calculation  is 
usually  smaller  than  that  of  a  Cl  calculation,  we  expect  to  obtain  potential 
energy  curves  of  a  higher  quality.  The  total  energies  obtained  will  generally 


To  assess  the  accuracy  of  the  computed  cross  sections,  we  compare  our 

results  with  experimental  and  other  theoretical  cross  sections  for  the  X^E^-to- 

B*£+  and  X-to-B' ^£+  transitions  in  Ho.  Cross  sections  for  transitions  to 
u  u  * 

bound  vibrational  states  are  calculated  in  addition  to  those  for  dissociative 
transitions.  The  experimental  and  previous  theoretical  results  do  not  include 
treatment  of  vibrational  and  rotational  effects.  Therefore,  we  compare  cross 
sections  in  which  the  rotational  degrees  of  freedom  are  treated  as  degenerate 
and  excitation  is  from  the  ground  vibrational  state.  In  general  the  agreement 
of  our  cross  sections  with  those  from  previous  theoretical  calculations  is 
very  good. 

The  experimental  cross  sections  for  the  X-to-B  transition  to  bound  states 
are  smaller  than  the  theoretical  ones  by  about  a  factor  of  2.  Most  of  this 
discrepancy  can  be  attributed  to  problems  with  resolving  the  measured  signal 
at  small  scattering  angles. 

The  theoretical  cross  sections  for  direct  dissociation  to  produce  H(2s) 
via  the  X-to-B'  transition  are  smaller  than  the  experimental  cross  sections 
by  factors  of  3  to  4.  However,  the  experiment  measures  the  total  H(2s)  formed; 
it  can  not  distinguish  products  of  direct  dissociation  through  the  B'  state 
from  products  of  excitation  through  other  electronic  states  or  indirect  proc¬ 
esses,  such  as  cascading  from  higher  electronic  states.  The  agreement  of 
the  IP  method  with  other  theoretical  methods  is  very  encouraging,  and  the 
comparisions  with  experiment  are  consistent  with  our  understanding  of  the 


processes . 


This  depends  strongly  upon  the  initial  and  final  state  potential  energy  curves 
and  the  transition  dipole  moment  between  the  two  states  but  only  weakly  upon 
the  interaction  of  the  incident  electron  with  the  target  molecule.  The 
dynamical  threshold  energies  are  probably  accurate  to  within  0. 1-0.5  eV  for 
all  but  the  HgBr  system.  The  energetic  and  dynamical  threshold  energies  are 
shown  in  Table  2.  The  dynamical  threshold  will  be  most  accurate  for  systems 
in  which  the  energetic  and  dynamical  thresholds  are  close,  such  as  the  X  to 
B  transition  in  O2  and  S2.  The  behavior  of  the  IP  and  BA  cross  sections  for 
energies  very  near  threshold  is  not  accurate  for  electronic  transitions.  These 
cross  sections  can  be  in  error  by  an  order  of  magnitude  or  more  because  the 
approximations  in  the  theories  break  down,  exchange  interactions  become  more 
important,  and  important  resonance  effects  are  neglected. 

Although  the  absolute  magnitudes  of  the  integral  cross  sections  may  be 
in  error  by  as  much  as  an  order  of  magnitude  for  the  lowest  energies, 
predictions  of  the  dependence  of  the  cross  sections  upon  initial  vibrational 
and  rotational  states  should  be  reliabl;.  The  change  in  cross  sections  as 
a  function  of  the  initial  vibrational  and  rotational  state  is  controlled  by 
the  change  in  overlap  between  the  initial  and  final  wavefunctions  for  nuclear 
motion.  Thus  it  is  the  structural  factor  determined  by  the  potential  energy 
curves  and  dipole  transition  matrix  elements  which  determines  the  vibrational 
and  rotational  dependence  of  the  cross  sections  and  not  the  dynamics  of  the 

interaction  between  the  incident  electron  and  the  target  molecule.  Also  the 

\ 

change  in  the  dynamical  threshold  energies  with  initial  vibrational  state 
should  be  accurate  to  within  0. 1-0.5  eV. 

More  accurate  dynamical  methods  are  required  to  obtain  accurate  absolute 
cross  sections  near  threshold.  The  first-order  perturbation  treatment  of 
the  collision  actually  becomes  valid  at  very  low  translational  energies  for 
nonzero  impact  parameters.  To  see  this  consider  the  effective  potential  felt 
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by  the  Incident  electron  at  long  range.  This  potential  is  the  sum  of  the 
(positive)  centrifugal  term  and  a  multipole  expansion  whose  leading  term  is 
a  (negative)  quadrupole.  The  potential  is  negative  at  small  distances  then 
increases  to  a  maximum.  The  value  of  the  potential  at  the  maximum  and  the 
distance  increase  with  increasing  impact  parameter.  If  the  electron  energy 
is  small  enough  the  centrifugal  barrier  will  prohibit  the  electron  from 
penetrating  to  close  distances  and  the  first  order  perturbation  approximation 
to  the  interaction  will  be  valid.  For  electronic  excitations  this  criterion 
can  only  be  met  in  the  exit  channel  and  this  approach  will  not  be  useful. 
However,  the  BA  or  IP  methods  could  be  used  to  accurately  predict  the  threshold 
behavior  of  the  cross  sections  for  rotational  and  vibrational  excitation  cross 
sections  which  are  dominated  by  scattering  at  nonzero  impact  parameter. 

The  only  methods  capable  of  predicting  accurate  cross  sections  at  threshold 
energies  are  those  based  upon  an  accurate  quantum  mechanical  treatment  of 

the  scattering  system.  These  include  close-coupling  calculations  and  the 
distorted  wave  Born  method  which  is  very  similar  to  a  two-Btate  close-coupling 
calculation.  The  major  limitation  in  the  distorted  wave  method  1b  the  improper 
treatment  of  resonance  behavior  which  iB  very  important  at  low  energies.  The 
limitations  in  the  close-coupling  calculations  are  associated  with  numerical 
difficulties  in  the  calculations.  The  major  difficulties  are  the  accurate 

inclusion  of  exchange  interactions,  the  related  difficulty  of  including 
electron-capture  states  in  the  expansion  of  the  total  wavefunction,  and  the 
convergence  of  the  cross  sections  with  respect  to  the  number  of  channels. 

Local  approximations  to  the  electron  exchange  terms  have  been  developed  and 
have  been  shown  to  be  extremely  easy  to  use  and  also  very  accurate  for  energies 
above  10  to  20  eV.  These  methods  have  not  been  adequately  tested  for  low 
energies  but  may  be  accurate  to  within  about  a  factor  of  two  in  the  cross 

section  for  threshold  energies.  This  could  be  easily  tested  for  some  atomic 


cases  for  which  accurate  quantum  scattering  calculations  have  been  performed 
down  to  very  low  energies.  Once  a  local  approximation  to  the  exchange  is 
included,  electron-capture  states  can  be  routinely  added  to  the  basis  without 
further  complications  and  resonance  features  can  be  accurately  described. 
Cross  sections  are  known  to  converge  slowly  with  increasing  basis  set  size. 
Rather  than  using  large  basis  sets,  accurate  calculations  can  be  obtained 
using  small  basis  sets  in  which  the  electronic  wavefunctions  for  the  target 
molecule  are  allowed  to  distort  as  the  incident  electron  approaches.  These 
polarization  effects  have  been  effectively  included  in  calculations  of  elastic, 
and  vibrationally  and  rotationally  inelastic  cross  sections,  but  have  not 
been  used  in  calculations  of  electronic  excitations. 

In  summary,  the  IP  cross  sections  for  electron-impact  dissociation  for 
energies  a  few  eV  above  threshold  are  probably  accurate  to  within  factors 
of  2  to  5  for  most  of  the  transitions  studied  here.  The  major  exception  is 
the  X-to-II  transition  in  S2  which  may  be  in  error  by  an  order  of  magnitude 
or  more.  The  cross  sections  are  qualitative  at  best  for  energies  near 
threshold.  For  these  energies  resonances  can  drastically  alter  the  shape  of 
cross  sections  and  they  are  not  included  in  the  IP  method.  The  effect  of 
increasing  the  initial  vibrational  state  should  be  well  predicted  and  the 
threshold  energies  and  their  dependence  upon  initial  vibrational  state  should 
be  accurate  to  within  0.1  to  0.5  eV. 


•  •  1 
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Comparison  with  Other  Dissociation  Processes 


A.  H2 


Dissociative  attachment  (DA)  is  the  predominant  dissociation  mechanism  in 
H2  below  5  eV.  The  process  is  understood  to  involve  resonant  excitation  of  low 
lying  electronic  states  of  H2,  which  subsequently  dissociate  into  H  and  H“.  A 
strong  dependence  of  the  cross  section  on  molecular  vibration  has  been  theoretically 
predicted-*  and  experimentally  verified^*^.  Theoretical  prediction  of  a  large 
rotational  enhancement^  was  not  observed  in  experiments^. 

The  DA  process  has  also  been  studied  experimentally  in  both  H2  and  D2  by 
Allan  and  Wong®,  who  investigated  with  refined  techniques  the  effect  of  both  mole¬ 
cular  vibrational  and  rotational  energy  on  the  dissociation  cross  section.  They 
concluded  that  there  is  a  strong  vibrational  effect  in  which  the  threshold  cross 
section  increases  by  four  orders  of  magnitude  for  vibrational  excitation  from 
v  ■  0  to  v  *  4.  They  determined  that  the  rotational  energy  dependence  is  much 
smaller,  and  is  described  well  by  the  resonance  model  calculations  of  Wadehra 
and  Bardsley^"** .  The  cross  section  increases  by  approximately  a  factor  of  five 
from  j  *  0  to  j  *  7.  Thus  as  a  function  of  total  internal  energy,  Initial  vibration 
is  much  more  effective  than  rotation  in  increasing  the  efficiency  of  dissociation 
by  this  mechanism.  However,  a  more  detailed  analysis**  indicates  that  it  is  the 
total  amount  of  available  internal  energy  that  determines  the  dissociation  rate, 
and  that  rotational  energy  is  efficient  provided  the  molecules  are  sufficiently 
rotationally  hot. 

The  resonance  dominated  dissociative  attachment  process,  important  below 
5  eV,  rapidly  diminishes  above  that  energy  to  values  below  10“2*  cm^.  Above  15 
eV,  a  major  channel  is  direct  dissociation  to  neutral  species.  This  process  is 
the  thrust  of  the  research  discussed  in  this  report.  As  for  DA,  the  direct  dissoci- 


ation  process  is  enhanced  by  vibrational  energy,  but  to  a  lesser  extent.  A  factor 


of  four  increase  in  the  cross  section  is  found  in  going  from  v  =  0  to  v  =  4.  In 
contrast  to  DA,  there  is  a  negligable  effect  due  to  rotational  energy.  The 
threshold  energies  decrease  by  the  vibrational  state  energy  change  in  going  to 
higher  V.  Since  we  are  already  at  fairly  high  energies  (above  20  eV)  these  incre¬ 
ments  represent  small  changes  in  total  energy. 

B.  02 

The  effect  of  initial  vibrational  energy  upon  the  dissociation  of  02  by 

electron  attachment  has  been  studied  both  experimentally^*^  and  theoretically^. 

In  the  following  figure  the  DA  cross  sections  are  shown  for  thermal  distributions 

of  vibrational  states  at  300  and  1400  K.  The  cross  sections  shown  are  those 

calculated  by  O'Malley,  which  agree  very  well  with  experiment.  At  300  K  over 

997.  of  the  02  molecules  are  in  their  ground  vibrational  state  and  at  1400  K 

approximately  807.  are  in  vj  s  0  and  16£  are  in  vj  “  1.  Therefore,  the  300  K  DA 

cross  sections  are  dominated  by  dissociation  from  vj  =  0,  whereas  the  DA  cross 

sections  for  1400  K  are  the  averages  of  cross  section  for  a  few  vibrational  states. 

Vibrational  excitation  of  the  02  molecule  gives  a  large  enhancement  of  the  cross 

section  at  threshold  energies.  The  cross  section  curves  are  shifted  to  lower 
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energies  by  about  1  eV  as  is  increased  from  0  to  1.  This  is  reflected  in 
a  shift  of  the  cross  section  curve  at  1400  K  to  lower  energies  relative  to  300 
K. 

For  comparison,  impact  parameter  cross  sections  for  direct  dissociation  through 
the  B  state  are  also  shown  for  ■  0,  1,  and  2.  The  IP  cross  sections  are  shifted 
to  lower  energies  by  an  energy  approximately  equal  to  the  vibrational  energy  spacing 
in  02  (about  0.2  eV)  as  is  increased  from  0  to  1.  The  different  behavior  in 
these  cross  sections  compared  to  the  DA  ones  at  threshold  arises  from  a  difference 
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in  the  mechanisms  for  the  two  processes.  The  nonionizing  dissociation  process 
occurs  by  excitation  from  the  vibrational  level  of  the  ground  electronic  state 
to  the  continuum  of  the  B  state  with  an  energetic  threshold  of  about  6.7  eV  for 
®  0.  The  IP  cross  sections  at  threshold  energies  are  controlled  by  the  overlap 

between  the  continuum  wavefunction  and  the  initial  vibrational  wavefunction.  In 

this  case,  the  overlap  does  not  depend  strongly  upon  the  initial  vibrational  energy 
level,  there  being  only  a  slight  enhancement  of  the  V£  =  l  and  2  cross  sections 
over  those  for  V£  ■  0  at  threshold  energies.  Therefore,  for  the  nonionizing  process 
the  thresholds  decrease  by  approximately  the  vibrational  energy  level  spacing 
as  vj  increases.  The  dissociative  attachment  of  0£  is  postulated  to  occur  through 
the  2nu  state  of  0^  with  a  threshold  of  3.6  eV  for  *  0.  The  cross  section 
for  this  process  is  believed  to  be  quite  high'(~10**^  cm2)  and  the  magnitude  of 
the  DA  cross  section  is  determined  by  the  competition  between  autoionization  to 
form  neutral  O2  and  dissociation  to  form  0+0*.  The  major  effect  in  this  process 
is  the  change  in  the  probability  of  survival  in  the  ionic  state  as  the  vibrational 
state  is  changed.  Because  this  enters  through  an  exponential  factor  in  the  cross 
section  it  is  seen  to  have  a  much  more  drastic  effect  upon  the  DA  cross  sections 

than  a  simple  shift  to  lower  energies  by  the  vibrational  energy  level  spacing. 

O'Malley  has  shown  that  the  effect  of  rotational  excitation  is  minor  compared 
to  the  vibrational  effect  for  the  DA  cross  section*2.  This  agrees  with  the  conclu¬ 
sion  we  have  drawn  for  the  direct  dissociation  process. 

We  note  that  the  effect  of  including  exchange  effects  will  change  the  shapes 
of  the  cross  section  for  the  direct  dissociation  process  at  threshold  energies. 
However,  the  mechanism  of  direct  excitation  will  be  dominated  by  the  overlaps 
of  the  initital  and  final  wavefunctions  and  the  qualitative  trends  seen  in  the 
IP  calculations  are  still  expected.  The  effect  of  resonances  near  threshold  are 
more  complicated  and  can  alter  even  the  qualitative  nature  of  the  predictions. 
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HC1 


The  dependence  of  the  dissociative  attachment  cross  section  for  HCl  has  been 
studied  by  Allan  and  Wong^.  The  mechanism  for  this  process  is  similar  to  that 
discussed  for  O2.  Electron  attachment  is  postulated  to  occur  to  form  the  lowest 
h+  state  of  HC1~  which  has  a  threshold  for  dissociation  only  0.6  eV  above  the 
ground  vibrational  state  of  the  ground  electronic  state  of  neutral  HCl.  The 
vibrational  energy  level  spacing  of  HCl  is  about  0.35  eV;  therefore,  by  *  2 
the  DA  cross  sections  are  nonzero  down  to  zero  translational  energy.  A  dramatic 
enhancement  of  the  DA  cross  section  is  seen  as  the  temperature  is  raised  from 
300  to  1000  K.  For  1000  K  the  DA  cross  sections  extend  down  to  zero  energy  even 
though  the  relative  population  of  vj  =  2  to  “  0  is  0.037..  As  in  the  case  of 
dissociative  attachment  for  O2,  for  HCl  the  cross  sections  for  electron  attachment 
to  HCl  are  very  high  at  low  energies  and  roughly  independent  of  initial  vibrational 
state.  The  DA  cross  sections  are  then  determined  by  the  competition  between 
dissociation  and  autoionization  to  form  neutral  BpecieB.  The  survival  factor 
(probability  of  not  autoionizing)  is  strongly  dependent  upon  the  vibrational  state 
and  enters  the  DA  cross  section  through  an  exponential  factor. 

The  structure  seen  in  the  DA  cross  sections  for  HCl  are  not  caused  by  the 
DA  process  itself  but  by  the  competition  of  dissociative  attachment  with  vibrational 
excitation  of  HCl.  At  energies  at  which  excited  vibrational  states  of  neutral 
HCl  become  accessible,  the  vibrational  excitation  cross  sections  become  very  large 
and  caused  a  dip  in  the  DA  cross  sections. 

The  impact  parameter  cross  section  for  dissociation  of  HCl  through  the  A 
state  are  shown  in  Appendix  E.  Near  threshold  the  effect  of  increasing  the  initial 
vibrational  level  is  quite  dramatic.  This  is  because  the  overlap  of  the  initial 
and  final  wavefunctions  for  energies  near  threshold  is  very  sensitive  to  v^.  Also, 
the  structure  seen  in  the  cross  section  near  threshold  is  a  reflection  of  the 
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structure  in  this  overlap  caused  by  the  nodes  in  the  excited  state  wavefunctions 
for  the  initial  electronic  state.  The  origin  of  this  structure  is  very  different 
from  that  for  the  structure  in  the  DA  cross  sections. 

As  mentioned  in  the  previous  section,  the  IP  method  is  not  valid  near  threshold 
but  the  dependence  of  the  cross  sections  upon  vibrational  energy  level  should 
be  reliably  predicted.  However,  near  threshold  resonance  effects  can  be  important 
and  can  significantly  alter  these  predictions. 


7. 


Recommendations 


The  cross  sections  presented  in  this  report  are  suitable  for  kinetic  modelling 

of  electron  dissociation  processes  in  plasmas  at  energies  above  threshold.  Since 

electron* impact  excitation  and  dissociation  of  molecules  plays  an  important  role 

in  laser  plasmas,  plasma  switches,  and  general  gas-phase  electric  discharges, 

it  is  important  to  develop  theoretical  methods  that  are  reliable  near  threshold. 

We  have  examined  possible  approaches,  and  conclude  that  a  more  rigorous  treatment 

using  discrete  variable  methods  is  promising  and  should  be  pursued. 

Extensions  of  the  impact  parameter  method  to  treat  polyatomic  molecules  is 

possible  by  using  wavepacket  methods  or  the  eikonal  approximation.  Both  methods 

14  15 

have  been  successfully  applied  to  photodissociation.  ’  Since  selective 

vibrational  excitation  of  the  parent  polyatomic  molecule  can  lead  to  changes  in 
the  ratio  of  product  fragments,  this  process  should  be  investigated  as  a  possible 
means  to  rapidly  switch  a  gas  discharge. 
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APPENDIX  A 


An  improved  impact  parameter  method  for  electronic  excitation  and 
dissociation  of  diatomic  molecules  by  electron  impact 


Michael  J.  Redmon,  Bruce  C.  Garrett, 
Lynn  T.  Redmon,  and  C.  W.  McCurdy* 


Chemical  Dynamics  Corporation 
1550  West  Henderson  Road 
Columbus,  OH  43220 


Permanent  address:  Department  of  Chemistry,  Ohio  State  University, 

Columbus,  Ohio  43210 


The  impact  parameter  method  for  electron  impact  excitation  of  diatomic 

ecules  is  reformulated  to  explicitly  treat  molecular  vibration  and  rotation 

to  permit  the  study  of  molecular  dissociation.  Applications  are  made 

optically  allowed  transitions  involving  the  X  ^£  +  .  B  and  B'  states 

g  u  u 

H2.  The  resulting  cross  sections  are  compared  to  other  theoretical 
culations  and  to  experimental  data.  This  method  is  applicable  to  heavy 
tomic  molecules  and  is  expected  to  be  useful  in  studying  trends  in 
ctronic  excitation  and  dissociation  cross  sections  associated  with 
iations  in  internal  energy. 
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X  £  (R)  — CkR(j  (kR)  -  tan(n )n  (kR))  (35) 
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AB  a  £ 


(36) 
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and  n  are  regular  and  irregular  Ricatti-Bessel  functions,  respectively, 
f 

is  the  phase  shift,  and  C  is  a  constant  obtained  by  normalizing  the  wave- 
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iction  to  unit  density  per  unit  energy 
>)]  is  replaced  by 


The  transition  energy  AE  (Eq. 
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d  the  final  expression  for  the  averaged  cross  section  is 
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ere,  in  analogy  to  Eq.  (28), 
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this  case,  the  final  rotational  angular  momentum  quantum  number  is 
terpreted  as  the  orbital  angular  momentum  quantum  number  If  of  the  separating 
oms .  We  refer  to  the  cross  section  defined  by  Eq.  (38)  as  the  1PVRD  cross 
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Invoking  the  Franck-Condon  approximation  in  Eq.  (12)  and  assuming  the 
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vibrational  states  are  degenerate  we  obtain 
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where  we  assume  the  sum  over  the  Franck-Condon  factors  is  unity  (this  is 

true  only  if  continuum  states  are  also  included  in  the  "sum").  This  expression 

thus  represents  excitation  to  both  bound  and  dissociative  states.  We  call 

the  cross  section  denoted  by  Eq.  (29)  the  IPVR  cross  section,  the  one  given 

by  Eq.  (31)  the  IPV  cross  section,  and  the  one  given  by  Eq.  (32)  the  IP  cross 
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section.  The  last  is  the  form  obtained  by  Hazi. 


B.  Bound- to-Continuum  Transitions 


Next  we  consider  the  modification  necessary  for  dissociative  transi¬ 
tions.  In  this  case  we  no  longer  have  a  vibrational  quantum  number  for  the 
final  state;  instead  we  have  the  continuum  energy  e.  Thus  w1?  replace  all 
occurrences  of  the  quantum  number  Vf  by  e  in  the  preceeding  equations.  The 
wave funct ions  for  the  continuum  states  are  solutions  to  the  equation 
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subject  to  scattering  boundary  conditions,  i.e., 
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fhere  ( T ro j- )  is  the  rotational  partition  function  for  vibronic  state  a^vf 
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kB  i 8  Boltzmann's  constant,  and  is  the  vibrational  quantum  number  of 

the  highest  bound  state  of  the  final  electronic  state.  The  nuclear  spin 

degeneracy  factor  dj^  is  unity  except  for  certain  classes  of  homonuclear 

diatomics  (see  Table  1).  j™aX  is  selected  by  practical  consideration. 

The  sum  over  final  rotational  states  has  nonzero  contributions  only 

for  jf  equal  to  jfi  or  jj+1,  (jf20).  In  the  low  temperature  limit,  only 

ataf 

the  jf=0  and  jf=l  states  contribute  to  o  (E,Trot).  Then  the  sums  over 

vi 

the  Clebsch- Gordan  coefficients  in  Eq.  (19)  can  be  evaluated  to  give  a  factor 
of  1/9  and  the  cross  section  Eq.  (29)  reduces  to 

max 


o  jC  f 

o  1  (E,T)  — ► 

Vi  Trot-0 


vf 

2 

V£  =  0 


2  2 
m  e 
e 

3ft4 


(2- 


6Af,0> 


IK 


iafai] 

vfvi 


3af“i 

Vfvi 


(31) 


I 


1 

i 

A 


44 


where 
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Y  = 

2 9.0  AE/mg(u^  +  u|) 

(26) 

E  is  the  relative  translational  energy  of  the  incident  electron,  m  is  the 
electron  mass,  l0  is  the  minimum  value  of  the  incident  electronic  orbital 
angular  momentum  quantum  number,  U£  (uf)  is  the  initial  (final)  electron 


speed,  and  the  transition  energy  AE  is  defined  by 


(27) 


The  dependence  of  the  dynamical  factor  on  molecular  vibrational  and  rotational 
quantum  numbers  is  only  through  the  transition  energy. 

The  cross  section  is  expressed  as  the  product  of  the  structural  factor 
for  the  molecule,  given  by  Eq.  (19),  and  the  dynamical  factor  for  the  electron, 
given  by  Eq.  (21)* 
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A.  Bound- to -Bound  Transitions 


The  cross  section  of  interest  is  obtained  by  averaging  Eq.  (28)  over 
a  Boltzmann  distribution  of  initial  rotational  states  characterized  by  a 
rotational  temperature  Trot  and  summing  over  all  final  vibrational  and 
rotational  states: 
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40 

inequality  .  The  structural  factor  is  zero  unless  j£  =  j£±l.  It  also  vanishes 
when  |  AA|  >1,  and  so  we  can  have  E  4— ►  E,  E  <—*•  n,  and  n  4 — ►  II  transitions  but 
not  E  4— *  A.  Extending  the  present  theory  to  permit  the  description  of  higher 
order  processes  would  require  including  additional  terms  in  the  multipole 
expansion  of  the  electron-target  interaction.  If  we  assume  degenerate 
rotational  states,  the  structural  factor  reduces  to 


c®  fai  = 

SV£Vi 


(2  -  6  )(2  -  6.  .) 

Af,0  A^O 


(20) 


Evaluation  of  the  dynamical  factor  follows  the  derivation  of  Hazi  exactly. 

By  assuming  that  the  trajectories  for  the  incident  electron  are  straight 

lines  and  energy  exchange  with  the  target  occurs  at  the  distance  of  closest 

approach  and  requiring  conservation  of  energy  and  angular  momentum,  the 
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dynamical  factor  can  be  expressed  in  terms  of  modified  Bessel  functions 
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K£(y)  and  modified  Struve  functions  S£(y)>  The  derivation  is  found  in 
Hazi's  paper,30  and  we  present  only  the  final  result  integrated  over  impact 
parameter: 
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We  have  also  used  the  relation  defining  the  displacement  vector  in  terms 
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of  the  spherical  harmonics 
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30 

to  simplify  the  time  dependent  integral.  This  form  for  the  probability 

results  from  the  present  restriction  to  dipole  allowed  transitions.  Thus 

even  with  the  inclusion  of  vibrational  and  rotational  degrees  of  freedom, 

we  still  obtain  a  form  analogous  to  Hazi's,  i.e.,  a  factorization  of  the 

transition  probability  into  the  product  of  a  dynamical  factor  depending  on 

projectile  electron  coordinates  (classical  trajectory)  and  a  unitless 

30 

structural  factor  depending  only  upon  properties  of  the  target  molecule. 

The  structural  factor  is  written  as 
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where  m  is  the  electron  mass.  It  contains  Clebsch-Gordan  coefficients  that 
e 

impose  optical  selection  rules  due  to  the  requirements  of  the  triangle 
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where  we  have  used  the  property 


j  *  *  m  -A.  j, 

DJi,  (R)  ■=  (-1)  1  i  DJi  .  (R) 

”rai_ Ai 


Making  use  of  the  integral  theorem  for  the  product  of  three  rotation 
matrices^*,  we  obtain 
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We  have  now  reduced  the  rotational  factor  to  a  simple  expression  involving 
the  product  of  two  Clebsch-Gordan  coefficients.  The  time-dependent 
perturbation  matrix  element  now  becomes 
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Substituting  Eq.  (16)  into  Eq.  (1),  summing  over  final  angular  momenta 
(electronic  and  rotational),  and  averaging  over  initial  angular  momenta  results 
in  the  following  form  for  the  transition  probability 
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V  ■!  V%  '  .■'X  1  J  ".'3  •  '•  ■  I  ' 


This  operator  is  a  spherical  tensor  of  rank  1,  so  (from  the  Wigner-Eckart 
39 

theorem  )  we  have 
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where  the  electronic  transition  dipole  moment  is 
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and  AA  =  Aj-A^.  The  transition  matrix  element  can  now  be  written  as 
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where  the  electronic  transition  dipole  moment  is  averaged  over  the  initial 
and  final  vibrational  wavefunctions 
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For  the  remaining  integral  over  Euler  angles,  we  use  the  Wigner  D-function 
representation  of  the  symmetric  top  wavefunctions,  Eq.  (4), 
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To  evaluate  this  expression,  we  use  the  usual  multipole  expansion 
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where  r<(r>)  is  the  lesser  (greater)  of  r'  and  rq'  .  We  keep  only  the 

asymptotic  (r1  >  rq')  terms  of  this  expansion  because  the  impact  parameter 

treatment  is  valid  only  for  collisions  that  do  not  penetrate  the  electron 
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cloud  of  the  target.  We  are  presently  interested  only  in  optical 

(dipole-allowed)  transitions  (for  which  the  dominant  contribution  is  from  fL  m 
1).  Therefore,  we  approximate  the  electron-target  interaction  by 
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Since  the  electron  coordinates  of  the  electronic  vavefunctions  employed 
in  the  matrix  element  given  by  Eq.  (5)  are  expressed  in  BF  coordinates, 
Eq.  (7)  must  be  transformed  to  BF  coordinates  before  the  integration  can 
be  carried  out.  Using  the  Wigner  rotation  matrices^*  we  obtain 
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Substituting  this  expression  into  Eq.  (5),  we  note  that  the  part  of  the 
resulting  expression  involving  integration  over  the  target  electron  coordinates 
is  the  matrix  element  of  the  spherical  representation  of  the  dipole  operator. 
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The  vibrational  wavefunctions  for  electronic  state  a  are  orthonormal  eigen¬ 
functions  of  the  equation 

l-(<k2/2»'AB)  d2/dR2  +  Va(R)  +  -R2j<j+1>/2MABR2]  X^(R)  =  eavi  x“j(R)  (3) 

where M._  is  the  reduced  mass  of  the  diatomic  molecule.  ea .  is  the  vibrational 
AB  vj 

energy  level,  and  V^(R)  is  the  potential  energy  curve  for  electronic  state 
a . 

The  symmetric  top  functions  are  represented  in  terms  of  Wigner 
D- functions^ 

i  ^  (2 j+1t  i*  - 

NJ  (R)  J  — _)  D\(R)  (4) 

mA  \  8tt  2  /  mA 

where  the  projectile  electron  is  assumed  to  be  initially  moving  along  the 

/v 

spaced- fixed  (SF)  Z  axis,  and  where  R  represents  the  two  (nontrivial)  Euler 
angles  relating  the  BF  frame  to  SF  axes.  This  normalization  assumes 
integration  over  all  three  Euler  angles. 

In  our  formulation  of  the  semiclassical  impact  parameter  method,  the 
projectile  electron  velocity  is  assumed  to  be  high  enough  that  the  collision 
time  is  short  compared  to  the  periods  of  molecular  vibration  or  rotation. 
Then  the  adiabatic  nuclei  approximation  is  expected  to  be  valid,  and  the 
time-dependent  perturbation  matrix  element  required  in  Eq.  (1)  is  obtained 
by  taking  the  matrix  element  of  the  projectile-target  interaction  between 
initial  and  final  molecular  states  represented  in  the  form  of  Eq.  (2).  If 
r'(t)  represents  the  classical  trajectory  of  the  projectile  electron,  and 
t  q  the  coordinates  of  a  target  electron  (the  primes  indicate  SF  coordinates) , 
then  this  matrix  element  can  be  expressed  as 
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II.  THEORY 


Hazi's  implementation  of  the  impact  parameter  method  employs  the  "fixed 
nuclei"  approximation^,  averaging  the  semiclassical  transition  probability^ 

p£i  •  I  1  i  v£l(e)  i> 

A  £=±A  £  Af=±Af 

over  all  molecular  orientations.  In  this  expression  gj  is  the  degeneracy 

of  the  initial  electronic  state,  and  Af  are  the  usual  projections  of 

electronic  angular  momentum  on  the  body-fixed  axis,  AEf^  is  the  transition 

energy,  and  Vf^(t)  is  the  time-dependent  matrix  element  of  the  perturbation. 
We  present  an  extension  of  the  theory  to  allow  treatment  of  those  processes 
that  require  consideration  of  nuclear  motion. 
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We  work  within  the  framework  of  the  adiabatic  nuclei  approximation 
and  employ  a  Born-Oppenheimer  factorization  of  the  wave  function  for  the 
target  molecule 


(x;  R)  -  <l>  (x;  R)  X“  .  (R)  Nj  <R> 

avA  -  aA  v  j  mA 


(2) 


The  function  represents  an  approximate  solution  (for  electronic  state 
a)  to  the  n-electron  Schrodinger  equation  for  the  electrons  in  the  molecule 
moving  in  the  field  of  the  "fixed"  nuclei,  x^j  is  the  vibrational  wavefunction 
for  the  (v,j)  rovibrational  level  of  electronic  state  a,  and  N^,  is  a 

IRA 

symmetric-top  function  for  a  diatomic  molecule  with  electronic  angular  momentum 
38 

component  A.  In  our  notation,  x  collectively  represents  the  coordinates 
of  the  target  electrons,  while  R  and  R  represent  the  magnitude  and  angular 
position  of  the  relative  position  vector  R.  All  target  electronic  coordinates 
are  referred  to  the  body-fixed  (BF)  frame  and  the  polar  z-axis  is  taken  along 
the  line  joining  the  nuclei  (R). 


products  of  Clebsch-Gordon  coefficients.  By  performing  suitable  averages, 
the  original  cross  section  formula  of  Hazi  is  obtained. 

The  appendixes  provides  the  relevant  computational  details,  while  Section 
IV  presents  the  results  of  the  cross  section  calculations  for  several  levels 
of  approximation  within  the  IP  method.  As  a  test  of  the  method  for  electronic 
excitation  to  bound  vibrational  states,  we  apply  the  IP  method  to  the  X^E*  to 

and  X^Eg  to  B^E^  transitions  in  H2  and  compare  with  previous 
32  12  16  23  30 

experimental J  and  theoretical  ’  '  *  studies.  We  also  apply  it  to  the 
direct  dissociation  of  ^(^Eg)  through  the  B^E^  state  and  compare  the  results 
with  experimental^’^  and  theoretical^  t26  studies.  Section  V  provides 
a  summary  and  the  conclusions. 
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wave  theories.  Thus  it  is  complementary  to  these  approaches  and  to  the 
high-energy  Born-type  approximations.  In  the  semiclassical  impact-parameter 
(IP)  method  the  motion  of  the  electrons  (including  the  incident  electron) 
is  separated  from  nuclear  motion  by  the  Born-Oppenheimer  approximation.  This 
allows  the  molecular  electronic  problem  to  be  treated  as  accurately  as 
necessary,  independently  of  the  scattering.  Hazi's  theory  assumes  degenerate 
rotational  states  and  averages  the  transition  probability  over  molecular 
orientations.  His  formulation  results  in  cross-sections  that  satisfy 
reciprocity,  and  this  represents  a  refinement  over  the  original 
impact -parameter  method.  However,  the  Franck-Condon  approximation  is  utilized, 
and  internuclear  distances  are  treated  as  static.  Thus  vibrational  and 
rotational  motion  are  neglected.  Because  of  these  approximations  in  Hazi's 
theory,  the  resulting  cross  sections  do  not  distinguish  between  bound 
( rovibrational )  and  unbound  (continuum)  channels.  When  the  method  is  extended 
to  include  nuclear  motion,  this  differentiation  can  be  made. 

In  Section  II  we  present  extensions  of  the  impact  parameter  model  to 
include  the  effects  of  initial  vibrational  and  rotational  excitation  upon 
electron  impact  induced  excitation  and  dissociation.  We  are  interested  in 
nonresonant  optically  allowed  electronic  excitations  of  diatomic  molecules. 
These  involve  dipole-allowed  spin-conserving  transitions.  These  excitations 
can  be  to  bound  rovibrational  levels  of  the  final  electronic  state  or  to 
parts  of  the  potential  energy  curve  above  the  dissociative  limit  of  the  final 
electronic  state.  We  consider  only  processes  of  direct  dissociation  through 
excited  electronic  states  rather  than  the  resonance-enhanced  dissociative 
attachment  process.  When  the  diatomic  molecule  is  treated  as  a  symmetric 
top,  electronic  and  rotational  angular  momentum  are  coupled;  however,  the 
resulting  cross  section  retains  a  separability  of  structural  and  dynamical 
factors  analogous  to  that  of  Hazi's  original  method  and  involves  simple 
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I.  INTRODUCTION 


Processes  involving  electron-molecule  collisions  are  important  in 
gas-laser  systems,  atmospheric  physics  and  photochemistry.  Because  of 

their  importance  in  laser  plasmas,!  collision-induced  electronic  excitations 
and  dissociation  by  electron  impact  have  been  the  subject  of  recent 
experimental  interest;2  however,  there  has  been  less  progress  in  the 
theoretical  description  of  these  processes.  Although  electron-molecule 
scattering  theory  has  been  significantly  advanced  in  the  areas  of  elastic 
scattering  and  rotational  and  vibrational  excitation, ^"5  the  al>  initio 
calculation  of  electronic  excitation  of  molecules  by  electron  impact  is 
relatively  new.  Early  work  in  this  area  includes  Born  or  other  "plane-wave" 
approximation  calculations®-!^  that  are  reliable  only  at  relatively  high 
energies.  In  addition  a  two-state  close-coupling  method!®-!®  has  been  applied 
to  electronic  excitation  and  dissociation  of  a  number  of  states  of  H2  and 
to  the  excitation  of  the  a^rig  state  of  N£!^.  In  these  calculations  electron 
exchange  is  correctly  included,  but  the  treatment  of  nuclear  motion  utilizes 
the  Franck-Condon  approximation.  The  theory  can  treat  both  singlet-singlet 
and  s inglet- triplet  transitions.  The  distorted-wave  Born  approximation 
employing  the  T-matrix^®-^®  or  R-matrix22  methods  has  also  received 

attention.  Most  applications  of  the  distorted  wave  methods  have  been  to 
processes  involving  excited  electronic  states  of  H2,  although  applications 
to  N2^  and  F22^  have  also  been  made.  The  method  is  most  applicable  to 
spin- forbidden  transitions  involving  short-range  interactions. 

Hazi  has  extended  the  semiclassical  impact-parameter  method2®-2^  to 
electronic  excitation  of  diatomic  molecules-!®.  The  formulation  neglects 
exchange,  and  is  best  suited  for  the  treatment  cf  optically  allowed  transitions 
that  require  many  partial  waves  in  conventional  close-coupling  or  distorted 


section.  The  limits  of  integration  in  Eq.  (38)  are  rigorously  given  by  the 
continuum  threshold 


Emin  “  V  (R  _  00 ) 


(40) 


and  by  energy  conservation 
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The  expression  for  dissociation  from  selected  initial  vibrational  states 
in  which  rotations  are  treated  as  degenerate  is  generalized  from  Eq.  (31) 
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where  the  structural  factor  is  given  by 
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and  the  dynamical  factor  is  given  by  Eqs.  (21)  -  (26).  The  cross  section 

defined  by  Eq.  (42)  is  called  the  1PVD  cross  section. 
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III.  COMPUTATIONAL  DETAILS 


The  potential  curves  for  the  X  ^E+,  B  *E+  and  B*  ^E+  states  of  Ho  are 

g  u  u  z 

43  44  it 

obtained  by  fitting  accurate  ab^  initio  *  and  experimental 4,3  data  in  a 
manner  similar  to  that  of  Bla^s  and  Truhlar^.  Over  the  region  for  which 
the  potential  data  is  available  the  points  are  fitted  to  a  cubic  spline 
function.  For  small  R  values  the  potential  is  written  as 


V  =  A  R"1  exp(-BR) 


where  the  parameters  A  and  B  are  determined  so  that  the  function  fits  the  data 
points  at  the  two  smallest  R  values.  For  the  region  of  R  values  greater 
than  the  spline-fit  region  the  potential  is  fit  to  the  form 


V  -  V0  -  C6  R-6  -  C8  R-8 


where  VQ  is  the  experimental  dissociation  energy  and  the  parameters  C5  and 
Cg  are  determined  so  that  the  function  fits  the  data  points  at  the  two  largest 
R  values.  The  cubic  spline  fit  is  restricted  to  match  the  function  values 
and  their  derivatives  at  the  end  points  of  the  data.  Our  fits  to  the  three 
potential  curves  used  in  this  work  are  presented  in  Fig.  1. 

The  transition  dipole  matrix  elements  as  a  continuous  function  of  R 
are  obtained  by  4-point  (3rd  order)  Lagrange  interpolation^  of  the  accurate 
ab  initio  data.^»^®  The  fits  to  the  matrix  elements  are  shown  in  Fig.  2. 

The  vibrational  eigenvalues  are  obtained  by  solving  Eq.  (3)  using 

A  Q 

Cooley's  algorithm  .  Once  an  eigenvalue  is  found,  the  values  of  the 
unnormalized  wavefunction  Xy  •($  are  stored  on  a  grid  of  evenly  spaced  R  values. 


The  normalization  factor  is  determined  by  extended  Simpson's  rule 


integration,^  and  the  stored  values  are  normalized.  Unnormalized  continuum 

wavefunctions  for  a  fixed  energy  e  are  obtained  by  Numerov  integration  and 

stored  on  the  same  grid  of  R  values  used  for  the  vibrational  wavefunction. 

The  asymptotic  boundary  conditions  are  applied,  and  the  normalization  constant 

of  Eq.  (35)  is  obtained  using  the  method  described  by  Chung,  Lin,  and  Lee. 

The  structural  factor  Saf0tl  .  is  obtained  from  Eq.  (19),  where  the 

vf)fviJi 

summation  over  the  Clebsch-Gordan  coefficients  is  straightforward  and  the 
square  of  the  dipole  transition  matrix  element  is  obtained  from  Eq.  (12). 
The  integration  in  Eq.  (12)  is  performed  using  extended  Simpson's  rule  on 
the  grid  of  R  values  on  which  the  wavefunctions  are  saved.  In  the  case  of 
bound- to-continuum  transitions,  the  quantum  number  Vf  is  replaced  by  the 
energy  eigenvalue  e,  but  the  computational  procedure  is  identical. 

The  dynamical  factor  is  calculated  from  Eq.  (21).  Evaluation  of  all 
of  the  modified  Bessel  and  Struve  functions,  including  the  finite  integral 
of  the  Struve  function,  is  accomplished  by  Lagrange  interpolation  from  the 
tables  of  Abramowitz  and  Stegun.^2  The  order  of  interpolation  used  is  that 
recommended  in  the  tables. 

The  impact  parameter  method  requires  that  a  minimum  value  of  the  orbital 

angular  momentum  of  the  incident  electron  be  supplied  to  evaluate  Eq.  (21). 

30 

For  the  comparisons  reported  here,  we  have  used  Hazi's  value  (bQ=1.7  ag). 

In  subsequent  applications,  we  define  the  minimum  value  of  the  impact  parameter 

bg  by  requiring  the  cross  section  from  Eq.  (32)  to  agree  with  that  of  the 

30 

Born  approximation  at  a  high  energy  .  The  choice  of  bg  is  discussed  further 
in  the  following  paper. ^0 

For  dissociative  processes,  the  integral  over  the  final  energy  is 
performed  using  repeated  Gauss-Legendre  integration.  The  total  range  from 
emin  to  Gmax  need  not  be  considered;  only  a  limited  range  of  e  values 
contributes  significantly  to  the  cross  sections  for  all  incident  energies. 
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IV.  RESULTS 


In  this  section  we  report  the  results  of  calculations  of  various 
electron- impact  processes  in  using  the  theoretical  impact-parameter  methods 
of  Section  II.  The  calculations  permit  comparison  with  experimental  and 
other  theoretical  results.  We  have  studied  excitation  and  dissociation  of 
ground-state  via  the  B  and  B'  states  and  the  effect  of  varying  initial 
molecular  vibrational  and  rotational  energy  on  these  processes.  The  latter 
studies  are  the  first  of  their  kind  for  electron  impact  dissociation  to  neutral 
species. 

In  Table  II  we  present  cross  sections  for  electron  impact  excitation 
from  the  ground  (X)  state  to  the  B  state  of  H^Cv^O).  Our  IPV  results  are 
compared  with  other  theoretical  calculations  and  with  experiment.  The 
comparison  is  also  shown  in  Fig.  3.  Because  the  level  of  theory  in  our  IPV 
calculations  is  meant  to  be  comparable  to  the  treatment  of  vibration  in  the 
majority  of  other  theoretical  methods,  we  use  the  Franck-Condon  approximation 
for  the  comparisons  in  Table  II.  Our  most  accurate  method  employs  explicit 
numerical  integration  of  the  dipole  transition  moment  over  the  vibrational 
wavefunctions  of  the  target  molecule. 

The  agreement  of  the  Born  results  with  our  IPV  cross  sections  in  Table 
II  is  good  at  energies  above  50  eV  but  poorer  at  lower  energies,  where  the 
Born  approximat  ion  is  known  to  break  down.  The  fact  that  the  two  seta  of 
Born  results  appear  to  have  a  different  low  energy  behavior  is  somewhat 
disturbing.  The  difference  is  larger  between  the  two  Born  cross  sections 
than  that  found  in  comparing  the  smaller  Born  cross  section  with  the 
Born-Ochkur  result  (which  included  exchange).  Fliflet  and  McCoy^  attribute 
this  difference  to  the  use  of  different  final-state  wavefunctions.  The  IPV 
results  are  in  reasonable  agreement  with  the  distorted  wave  cross  sections, 
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being  somewhat  larger  at  the  lower  energies.  The  close-coupling  results 
are  also  smaller  than  the  IPV  cross  sections,  but  the  overall  agreement  of 

the  present  impact  parameter  method  with  both  of  these  more  exact  methods 

is  within  an  acceptible  15  percent  above  25  eV.  The  distorted  wave  and 
close-coupling  calculations  include  exchange,  which  is  not  treated  in  the 
impact-parameter  theory.  Exchange  effects  could  be  included  in  the  IP  method 
by  the  use  of  effective  exchange  potentials. 51  However,  one  should  then 
worry  about  higher  order  terms  in  the  expansion  Eq.  (6). 

The  experimental  results  of  Srivastava  and  Jensen^  lie  about  a  factor 
of  two  below  all  of  the  theoretical  calculations.  Their  experiment  did  not 
directly  measure  the  integrated  cross  section  to  all  final  vibrational  states; 
it  only  detected  radiation  from  the  v'=2  state.  The  quantity  they  measured 
was  the  differential  cross  section,  which  agreed  well  with  the  distorted 

wave  calculations  of  Fliflet  and  McCoy  except  in  the  low  angle  region,  where 
the  experimental  cross  section  was  about  a  factor  of  two  smaller  than  the 
calculations.  The  cross  section  summed  over  all  final  vibrational  states 
was  estimated  using  the  Franck-Condon  principle.  The  differences  between 
the  theoretical  and  experimental  cross  sections  in  Table  II  probably  stem 
from  the  smaller  values  for  low  angle  scattering  obtained  in  the  experiment. 
The  good  agreement  above  threshold  among  all  of  the  theoretical  methods 

suggests  that  the  predicted  values  of  the  excitation  cross  sections  in  this 
energy  region  for  the  X-to-B  transition  in  are  reliable. 

In  Table  III  we  compare  the  results  for  this  transition  of  various 
treatments  of  molecular  vibration  within  the  impact  parameter  method.  The 
IP  method  of  Hazi,^0  which  assumes  vibrational  degeneracy,  underestimates 
the  IPV  cross  sections  obtained  with  an  accurate  treatment  of  vibration  (shown 
in  the  right-most  column)  by  approximately  ten  percent.  The  use  of  IPV  with 
the  Franck-Condon  approximation  shows  the  improvement  over  IP  that  can  be 
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achieved  by  including  vibrational  motion.  The  IP  method  implicitly  assumes 
a  sum  over  bound  and  continuum  vibrational  levels,  while  the  IPV  method  only 
includes  bound  levels.  For  this  comparison,  however,  the  major  differences 
arise  from  errors  in  treating  these  bound  states  since  the  dissociation  cross 
section  is  small.  This  comment  applies  only  for  this  case  (v=Q).  As  will 
be  seen  below,  there  is  a  considerable  enhancement  of  the  dissociation  cross 
section  with  increasing  target  vibrational  excitation. 

Table  IV  presents  results  from  excitation  to  the  B'  state  of  H^.  For 
this  transition  there  are  fewer  calculations  for  comparison  and  (to  our 
knowledge)  no  experimental  results.  The  most  accurate  theoretical  calculations 
available  for  comparison  are  the  distorted  wave  results  of  Mu-Tao,  Lucchese, 
and  McCoy^.  These  results  include  bound  vibrational  states  through  a 
Franck-Condon  factor  in  the  differential  cross  section  so  that  the  most 
meaningful  comparison  is  with  our  IPV(FC)  calculations  in  the  fourth  column. 
The  agreement  with  the  present  IPV  cross  sections  is  typically  to  within 
thirty  to  forty  percent  of  the  distorted-wave  reBult.  The  agreement  between 
the  present  results  and  the  Born  calculations  of  Arrighini  et.  al.  is  poor, 
the  difference  being  a  factor  of  four  at  low  energies.  The  simple  IP  method 
overstimates  the  IPV  and  DW  cross  sections  by  approximately  a  factor  of  two. 


This 

is 

due  to 

the 

implicit 

sum  over  continuum  states  discussed 

above  for 

the 

IP 

method. 

The 

sum  over 

all  bound 

vibrat ional 

levels  of  the 

B'  state 

of 

the 

Franck- 

Condon 

factors 

involving 

ground-state 

H  is  0.48, 

so  that 

approximately  fifty  percent  of  the  IP  cross  section  comes  from  the  continuum. 
This  comparison  of  the  IP  and  IPV  cross  sections  indicates  that  a  realistic 
treatment  of  molecular  vibration  is  important  for  this  transition  due  to 
the  importance  of  dissociative  processes. 

In  Table  V  we  present  cross  sections  for  the  production  of  H(ls)  +  H(2s) 
via  excitation  to  continuum  levels  of  the  B'  state.  The  agreement  between 


the  various  theoretical  results  Is  good;  the  largest  difference  is  less  than 
twenty  percent  at  energies  above  the  maximum  in  the  cross  section.  The 
theoretical  results  all  lie  approximately  fifty  percent  below  the  experimental 
results  of  Vroom  and  de  Heer.  The  latter  have  been  scaled  by  0.8  as  suggested 
by  Mumma  and  Zipf^4  to  correct  for  molecular  radiation  observed  in  the 
experiment.  Because  the  experimental  cross  sections  contain  contributions 
from  other  states  that  produce  H(ls)  +  H(2s),  they  are  expected  to  be  larger 
than  the  theoretical  calculations.  Figure  4  illustrates  these  results  and 
shows  the  good  agreement  of  the  various  theoretical  treatments. 

In  Fig.  5  the  effect  of  increasing  the  energy  in  molecular  vibration 
(up  to  v=3)  for  the  bound-bound  X-B  transition  is  shown.  For  this  range 
of  vibrational  quanta  the  cross  section  doubles.  In  Fig.  6  the  corresponding 
vibrational  dependence  for  dissociative  excitation  is  seen  to  be  even  greater, 
although  the  cross  sections  are  smaller  than  for  nondissociative  excitation. 
The  enhancement  in  both  cases  results  form  the  larger  Franck-Condon  factors 
for  the  excited  vibrational  levels  of  the  X  state. 

In  Fig.  7  the  vibrational  dependence  of  the  cross  section  for 
nondissociative  X-B'  excitation  is  shown.  Instead  of  increasing  monotonically 
with  v,  the  v=l  and  v=2  cross  sections  reverse  order.  The  cross  sections 
for  dissociation  to  H(ls)  +  H(2s)  in  Fig.  8  are  monotonically  decreasing 
with  v,  in  contrast  to  all  of  the  other  processes  discussed  here.  Again, 
this  behavior  is  primarily  determined  by  the  magnitude  of  the  structural 
factor . 

The  impact  parameter  method  as  formulated  in  this  work  is  capable  of 
providing  information  on  final-state  vibrational  distributions.  In  Fig.  9 
we  show  vibrational  distributions  for  both  the  B  and  B'  states  resulting 
from  e lec t ron- impact  excitation  from  the  ground  vibrational  level  of  the 
X  state.  The  two  distributions  reflect  the  difference  in  the  number  of 
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vibrational  levels  supported  by  the  two  excited  electronic  states.  Both 
the  B  and  B'  distributions  have  maxima  for  v>0,  indicating  a  preference  for 
vibrationally  inelastic  processes. 

Table  VI  presents  the  rotational  temperature  dependence  of  the  X-B 
excitation  of  initially  in  the  ground  vibrational  state.  The  cross  section 
increases  by  approximately  forty-five  percent  over  the  temperature  range 
from  zero  to  20,000  K.  The  corresponding  dissociative  cross  section  (not 
shown)  to  form  H(ls)  +  H(2p)  increases  by  thirty  percent  over  the  same  range 

of  temperatures.  For  v=0  the  dissociation  cross  section  has  a  maximum  for 

-19  2  -17  2 

300  K  of  2.9x10  cm  at  50  eV  compared  to  5.6x10  cm  for  the  bound-bound 

excitation.  For  the  X-B  transition  bound-bound  transitions  dominate  over 

bound -continuum  processes  regardless  of  rotational  temperature  or  initial 

vibrational  temperature.  The  effect  of  initial  vibrational  energy  is  much 

more  significant  than  the  effect  of  rotational  temperature. 

In  Table  VII  the  nondissociative  excitation  ccross  section  for  the  X-B' 
transition  decreases  monotonically  with  temperature.  This  is  the  opposite 
of  the  behavior  seen  in  Table  VI  for  the  X-B  transition.  The  explanation 
lies  in  the  far  fewer  bound  rotational  states  that  can  be  supported  by  the 
B'  electronic  state.  At  higher  temperatures,  therefore,  the  tail  of  the 
Boltzmann  distribution  becomes  severly  truncated.  The  cross  sections  for 
dissociation  through  the  B  state,  shown  in  Table  VIII,  increase  monotonically 
with  temperature.  We  also  found  this  increase  with  temperature  for  the  X-B 
dissociation  cross  sections.  We  are  not  reporting  the  latter  since  they 
are  small  in  magnitude  for  v=0  (the  only  case  studied  when  rotational  motion 
was  Included)  and  the  trends  follow  those  for  the  B'  state. 


V.  CONCLUSIONS 


The  purpose  of  this  work  was  to  extend  the  impact  parameter  method  for 
diatomic  molecules  to  permit  a  realistic  treatment  of  vibrational  and 
rotational  motion  and  to  allow  the  study  of  dissociative  processes.  To  test 
the  theory,  applications  were  made  to  nondissociative  excitation  of  the  B 
state  of  H2  and  to  dissociation  to  H(ls)  +  H(2s)  via  the  B'  state.  Agreement 
with  previous  theoretical  results  is  generally  good.  Comparison  with 
experiment  was  hampered  by  the  lack  of  direct  experimental  information  for 
the  processes  we  studied.  However,  the  magnitudes  of  the  theoretical  cross 
sections  are  very  consistent  with  our  understanding  of  the  processes  probed 
by  the  experiments.  New  results  from  this  work  include  the  dependence  on 
initial  vibrational  state  and  rotational  temperature  of  nondissociative  and 
dissociative  processes  involving  the  X,  B,  and  B'  states  of  H2.  While  the 
dependence  of  the  cross  sections  on  initial  rotational  temperature  is  modest, 
we  find  that  increasing  initial  vibrational  energy  can  have  a  dramatic  effect 
on  the  cross  sections.  There  is  to  our  knowledge  no  experimental  data 
providing  information  on  the  effect  of  internal  energy  on  the  processes 
reported  here.  However,  the  trends  we  find  concerning  the  effect  of  internal 
energy  on  dissociation  are  similar  to  those  seen  in  dissociative  attachment 

4  u  52’53 
processes  in  H2. 

In  summary,  the  impact  parameter  theory  presented  here  should  be  useful 
in  studying  the  effect  of  internal  energy  in  excitation  and  dissociation 
processes  in  diatomic  molecules.  It  is  a  first  order  theory,  and  thus  is 
most  valid  at  high  energies,  i.e.  at  some  distance  above  threshold.  It  is 
also  best  applied  to  qualitative  studies  of  integral  cross  sections. 


Applications  of  the  methods  developed  here  to  heavier  diatomic  molecules 
are  given  in  the  following  paper^  and  in  forthcoming  publications . ^ ^ 
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Fi'TUf'  A- 4 


ENERGY  (eV) 


T 


T 


represent  the  Born-Rudge  results  of  Ref.  12  and  the  triangles  represent 
the  distorted-wave  results  of  Ref.  23.  The  experimental  results  of 
Ref.  30  (multiplied  by  0.8  according  to  the  recommendation  of  Ref.  31) 
are  shown  as  bullets.  The  experimental  results  include  contributions 
from  states  other  than  B1  ^Ey. 

5.  Dependence  of  the  cross  section  (summed  over  final  vibrational  states) 

for  excitation  of  the  B  1-Ey  state  on  initial  vibrational  quantum  number 
(in  units  of  10”^  cm^ ) .  The  solid  curve  represents  v=0,  the  dashed  curve 
v*l,  the  long  dash-short  dash  curve  v=2,  and  the  long  dash  curve  v=3. 

6.  Dependence  on  initial  vibrational  quantum  number  of  the  electron  impact 

dissociation  cross  section  to  produce  H(ls)  +  H(2p)  via  the  B  ^Ej  state 
of  H2  (in  units  of  10-^®  cm^).  The  solid  curve  represents  v=0,  the  dashed 
curve  v««l,  the  long  dash-short  dash  curve  v=2,  the  long  dash  curve  v=5, 
and  the  long  dash-double  short  dash  curve  v=10. 

7.  Dependence  of  the  cross  section  (summed  over  final  vibrational  states) 

for  excitation  of  the  B'  *Ej  state  on  initial  vibrational  quantum  number, 
(in  units  of  10“^-®  cm^).  The  solid  line  represents  v=0,  the  dashed  line 
v=l,  the  long  dash-short  dash  curve  v=2,  and  the  long  dash  curve  v=3. 

8.  Dependence  on  initial  vibrational  quantum  number  of  the  electron  impact 

dissociation  cross  section  to  produce  H(ls)  +  H(2s)  via  the  B'  ^Ey  state 
of  H2  (in  units  of  10"^®  cm^).  The  solid  curve  represents  v=0,  the  long 
dash-short  dash  v=l,  the  long  dash  curve  v=2,  and  the  short  dash  curve 
v  =  3 . 

9.  Final  vibrational  distributions  (in  units  of  cm^ )  of  the  B  ^eJ  and  B'  *Ey 

states  of  H2  resulting  from  electron  impact  excitation  of  the  X  ^E*  (v=0) 

state . 
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FIGURE  CAPTIONS 


1.  Potential  energy  curves  for  the  three  states  of  H2  studied  in  this  work. 

The  X  ^Eg  and  B  ^-Ey  potentials  are  fits  to  the  calculations  of  Kolos 
and  Wolniewicz  (Ref.  40  and  Ref.  41,  respectively),  while  the  B*  poten¬ 

tial  is  a  fit  to  the  data  of  Spindler  (Ref.  42). 

2.  Dipole  matrix  elements  used  in  the  calculations.  The  matrix  element 
for  the  X  ^Eg  to  B  ^Ey  transition  is  a  fit  to  the  calculation  of  Wolniewicz 
(Ref.  44)  while  that  for  the  X  ^E|  to  B'  ^E^  transition  is  a  fit  to 
the  calculation  of  Ford  (Ref.  45). 

3.  Electron-impact  excitation  cross  sections  (units  of  lO-^  cm2)  for  the 

X  (v=0)  to  B  ^Ej  transition  in  H2.  The  solid  curve  represents  the 

results  (summed  over  final  vibrational  states)  of  the  present 
impact-parameter  IPV  cross  section  (Eq.  (31)]  with  the  Franck-Condon 
approximation  used  to  evaluate  the  vibrational  integrals.  The  circles 
and  bullets  are  the  resultB  of  Born  approximation  calculations  (Ref. 
16  and  Ref.  20,  respectively).  The  triangles  are  from  the  L2  distorted 
wave  calculations  of  Ref.  20.  The  asterisks  are  the  Born-Ochkur 
calculations  of  Ref.  16,  and  the  squares  are  the  results  of  two-state 
close-coupling  calculations.  The  vertical  lines  represent  error  bars 
for  the  experimental  results  of  Ref.  29. 

4.  Electron-impact  dissociation  cross  sections  (in  units  of  lO”*®  cm2)  for 

production  of  H(ls)  +  H(2s)  via  excitation  of  the  B'  ^Ey  (v=0)  stabe  of 
H2-  The  solid  curve  is  the  result  of  applying  the  IPVD  method  of  Section 
II. B  with  accurate  numerical  evaluation  of  the  integrals  over  vibrational 
wavefunctions .  The  dashed  curve  represents  the  IPVD  result  when  the 

Franck-Condon  approximation  is  used  to  evaluate  the  integrals.  The  squares 


; 
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TABLE  A-viII .  Rotational  temperature  dependence  of  the  croas  section  for 

dissociation  of  through  the  B'  state. 


Temperature  (Kelvin) 


E(eV) 

0 

4.2 

300 

1000 

5000 

20000 

20 

2.28 

2.28 

2.31 

2.40 

2.77 

2.96 

30 

3.06 

3.06 

3.11 

3.24 

3,77 

4.04 

40 

3.46 

3.46 

3.51 

3.65 

4.21 

4.49 

50 

3.52 

3.52 

3.57 

3.70 

4.25 

4.53 

60 

3.45 

3.45 

3.49 

3.62 

4.14 

4.41 

80 

3.19 

3.19 

3.23 

3.34 

3.81 

4.05 

2.92 

2.92 

2.96 

3.06 

3.48 

3.69 

150 

2.38 

2.38 

2.42 

2.50 

2.84 

m 

200 

2.02 

2.02 

2.04 

2.12 

2.40 

2.54 

250 

1.76 

1.76 

1.78 

1.84 

2.08 

2.20 

1.56 

1.56 

1.58 

1.63 

1.85 

1.95 

_18  2 

Cross  sections  are  in  units  of  10  cm  . 
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TABLE  A-VII.  Rotational  temperature  dependence  of  the  cross  section  for 
electronic  excitation  to  the  B'  state  of  H ^ • 


Temperature  (Kelvin) 


E(eV) 

0 

4.2 

300 

1000 

5000 

20000 

20 

2.10 

2.10 

2.10 

2.14 

1.87 

1.65 

30 

2.90 

2.90 

2.89 

2.94 

2.57 

2.27 

3.15 

3.15 

3.14 

3.19 

2.78 

2.44 

50 

3.13 

3.13 

3.12 

3.17 

2.75 

2.42 

60 

3.03 

3.03 

3.01 

3.06 

2.65 

2.32 

80 

2.82 

2.82 

2.74 

2.78 

2.41 

2.12 

100 

2.50 

2.49 

2.52 

2.18 

1.92 

Cross  sections  from  Eq. 


(29)  are  in  units  of  10 


18 


2 

cm  . 
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table  A-VI.  Rotational  temperature  dependgnce  of  the  cross  section  for 
electronic  excitation  to  the  B  state  of  H^. 


E(eV) 

0 

4.2 

Temperature 

300 

(Kelvin) 

1000 

5000 

20000 

15 

3.87 

3.87 

3.89 

3.94 

4.31 

5.53 

20 

4.47 

4.46 

4.49 

4.57 

5.15 

6.83 

25 

5.25 

5.24 

5.27 

5.37 

5.97 

7.63 

30 

5.61 

5.60 

5.63 

5.72 

6.31 

7.93 

35 

5.72 

5.70 

5.74 

5.82 

6.39 

7.91 

40 

5.68 

5.67 

5.71 

5.79 

6.32 

7.75 

45 

5.58 

5.57 

5.60 

5.68 

6.19 

7.54 

50 

5.45 

5.44 

5.47 

5.54 

6.02 

7.29 

60 

5.14 

5.14 

5.15 

5.23 

5.66 

6.81 

80 

4.55 

4.55 

4.57 

4.63 

5.00 

5.95 

100 

4.07 

4.06 

4.08 

4.14 

4.45 

5.27 

Cross  sections  from  Eq.  (29)  are  in  units  of  10  cm  . 


a 
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table  A-v.  Electron  impact  dissociation  cross  sections  for  production 
of  H(ls)  +  H(2s)  via  excitation  to  the  B'  state  of  H^. 


E(eV) 

icsn:ttB:ssecxBCBB:i 

BRb 

DWC 

IBBBBttVBBetSBBtSBCI 

IPVD(FC)d 

BBBBBBEII 

IPVD® 

EXP, 

20 

1.98 

1.48 

2.32 

2.24 

30 

3.23 

3.12 

3.02 

3.91 

4.46 

3.54 

3.41 

50 

3.60 

3.47 

6.17 

60 

3.75 

4.41 

3.52 

3.40 

5.94 

70 

3.40 

3.28 

80 

3.41 

3.26 

3.14 

5.50 

90 

3.12 

3.00 

3.09 

3.41 

2.99 

2.88 

5.08 

150 

2.50 

2.44 

2.35 

4.19 

2.10 

2.07 

1.99 

3.62 

250 

1.80 

1.73 

3.10 

1.60 

1.54 

2.69 

a  -18  2 

Cross  sections  are  in  units  of  10  cm  .  Results  are  for  0  R. 

b  Born-Rudge  results  of  Ref.  12. 

c  2 

L  distorted -wave  results  of  Ref.  26. 

d  Impact  parameter  results  with  IPVD  cross  section  from  Eq.  (42) 
evaluated  using  the  Franck-Condon  approximation  for  vibrational  integrals. 

e  Impact  parameter  results  with  IPVD  cross  section  from  Eq.  (42) 
with  the  vibrational  integrals  evaluated  numerically. 

*  Experimental  results  of  Ref.  33  scaled  by  0.8  according  to 
the  recommendation  of  Ref.  34.  The  experimental  cross  section 
contains  contributions  from  states  other  than  B'. 
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table  A- IV.  Cross  sections  for  electron  impact  excitation 

from  the  X  to  the  B'  state  of  Hj. 


E(eV) 

b  0 

Born  DM 

IBBBBBCSB 

IP" 

ipv(fc)6 

IPVf 

20 

8.88  1.7 

4.03 

1.93 

2.14 

30 

3.9 

5.49 

2.66 

2.95 

40 

4.8 

6.19 

2.89 

3.21 

50 

9.41  4.9 

6.28 

2.88 

3.19 

60 

4.7 

6.14 

2.78 

3.08 

70 

4.4 

5.91 

2.66 

2.94 

80 

4.1 

5.79 

2.53 

2.80 

90 

3.8 

5.42 

2.41 

2.67 

100 

6.75  3.4 

5.18 

2.30 

2.54 

a  —18  2 

Cross  sections  are  in  units  of  10  cm  .  Results  are  for  0  K. 

k  Born  Results  of  Ref.  15. 

c  2  . 

L  distorted-wave  results  interpolated  from  Fig.  4  of  Ref.  26. 

d  Impact  parameter  IP  results  from  Eq.  (32)  corresponding 

to  the 

original  method  of  Hazi  (Ref.  30]  which 

assumes 

vibrational  degeneracy. 

Continuum  contributions,  implicit  in  this  method 

,  are  large 

for  this  case. 

0 

Impact  parameter  results  with  the  IPV  cross 

section  of 

Eq.  (31). 

The  vibrational  integrals  are  evaluated 

using  the  Franck-Condon  approximation 

^  Impact  parameter  results  with  the  IPV  cross 

section  of 

Eq.  (31). 

The  vibrational  integrals  are  evaluated 

using  numerical  integration. 
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table  a- hi.  A  comparison  of  various  impact-parameter  approximations 

r _ 1. 1 _ «  a _ _ '  l  • _ _  r  tt  ™ 


for  the  X  to  B  excitation  of  Hj. 


~1 7  2 

a  Cross  sections  are  in  units  of  10  cm  .  Results  are  for  0  K. 

k  ip  results  of  Eq.  (32)  corresponding  to  Hazi's  method  of  Ref.  30. 
Continuum  contributions  (small  in  this  case)  are  implicitly  included. 

c  IPV  results  of  Eq.  (31)  with  the  Franck-Condon  approximation  used 
to  evaluate  integrals  over  the  vibrational  wavefunctions. 

d  IPV  results  of  Eq.  (31)  with  the  integrals  over  vibrational 
wavefunctions  evaluated  explicitly. 


table  A- ii .  A  comparison  of  the  present  impact-parameter  results  for  the 

electron-impact  excitation  of  the  B-state  of  with  other  results. 


E(eV) 

ccb 

CBBCBBRI 

Born  IC 

B0d 

DW* 

Born  11^ 

IPV8 

e  h 
Exp. 

15 

1.05 

3.60 

3.61 

1 .4  ±  0.4 

20 

3.09 

5.23 

4.10 

1,9  ±0,6 

25 

4.31 

6.66 

5.31 

4.12 

5.64 

4.84 

30 

4.46 

5.69 

5.19 

2.0  ±  0.6 

40 

4.93 

5.43 

5.27 

2.3  ±  0.7 

50 

4.71 

5.55 

5.14 

4.84 

5.06 

5.06 

2.8  ±0.8 

60 

4.49 

4.70 

4.79 

1 .9  ±  0.6 

75 

4.09 

4.55 

4.36 

4.37 

100 

3.58 

3.87 

3.76 

3.80 

a  - 
Cross 

sections 

are  in 

units  of 

10  ^  cm 

2 

.  Results 

are  for 

0  K. 

Two-state  close-coupling  resits  of  Ref.  16. 
c 

Born  results  of  Ref.  16. 

d  Born-Ochkur  results  of  Ref.  16. 
e  2 

L  distorted- vave  results  of  Ref.  23. 

^  Born  results  of  Ref.  23. 

8  IPV  results  of  the  present  work  from  Eq.  (31)  with  Franck-Condon 
evaluation  of  vibrational  integrals. 

Experimental  data  from  Ref.  32. 
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TABLE  A- 1.  Nuclear  spin  degeneracy  factors  required  for  homonuclear 
diatomic  molecules  in  sigma  states. 


BBCSSBBBtiC&Bl 


EtBCttCBC 


d. 

3 


+ 

g  or  u 


and  21+ 1  even 


g 


+ 

or  u 


and  21+ 1  even 


g 


+ 

or  u 


or 

and  21+ 1  odd 


+ 

g  or  u 


or 

and  21+ 1  odd 


even  j  ( 1+1 )( 21+1 )  1(21+1) 

odd  j 


1(21+1) 


(1+1) (21+1) 


ENERGY  (eV) 


ENERGY  (eV) 


APPENDIX  B 


MCSCF  Calculations  of  the  Potential  Curves 
and  Dipole  Transition  Moments  for  Low-Lying  States 

of  O2  and 


by 

L.  T.  Redmon  and  R.  N.  Diffenderfer 


Chemical  Dynamics  Corporation 
1550  West  Henderson  Road 
Columbus,  Ohio  43220 


This  paper  describes  multiconfiguration  self-consistent  field  (MCSCF) 
calculations  of  potential  energy  curves  for  low-lying  states  of  O2  and  S2. 
The  states  of  interest  are  those  to  which  transitions  from  the  ground  state 
are  dipole-allowed.  The  dipole  transition  moments  are  also  presented  here. 
The  results  of  dynamical  calculations  of  electron  impact  cross  sections 
employing  the  present  data  are  given  elsewhere^-^. 

02 

The  potential  curves  for  O2  were  obtained  using  a  basis  set  of  44  atomic 

orbitals.  This  contracted  set  was  designed  by  starting  from  Dunning's  4s 

3p  basis  set^,  adding  d  orbitals  with  exponents  of  1.21  and  0.2,  and  adding 

a  set  of  diffuse  p  orbitals  (with  orbital  exponents  of  0.03)  to  describe 

3  3  - 

the  Rydberg  character  of  these  states.  The  first  17  and  X  states  have 

u  u 

been  shown^"^  to  become  Rydberg-like  at  short  internuclear  distances  as  a 

result  of  avoided  crossings  of  Rydberg  and  valence  type  diabatic  states.  In 

3 

the  case  of  the  11^  states  a  multiconfiguration  (MC)  space  consisting  of 

2  2  2  4  1  1 

the  48  terms  resulting  from  the  full  excitation  of  the  2o  2o  3o  lii  lir  3o  val- 

8  u  g  u  g  u 

ence  reference  configuration  into  the  Itt  and  3o  orbitals  is  sufficient  to  de- 

g  u 

scribe  the  purely  repulsive  potential.  The  MCSCF  method  is  even  able  to 
describe  this  state  in  its  Rydberg  region  because  of  its  ability  to  change 
the  orbitals  from  valence  to  Rydberg  type  in  order  to  satisfy  its  optimal 
variational  requirements.  However,  the  method  would  not  be  expected  to 
describe  both  the  upper  and  lower  states  in  the  avoided  crossing  region  unless 
the  MC  space  is  expanded. 

3 

Since  we  have  not  obtained  exact  (full  Cl)  results  for  the  X  and  IT 

u 

states,  the  orbitals  in  terms  of  which  they  are  described  are  not  perfectly 
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orthogonal  to  each  other.  Our  transition  dipole  moment  codes  require 
orthogonality.  Therefore,  there  are  two  ways  of  obtaining  matrix  elements 
between  these  states.  We  can  use  the  (optimal)  ground  state  orbitals  to 
obtain  a  nonoptimal  (though  variational)  wavefunction  for  the  pi  state  and 
then  evaluate  the  moment,  or  we  can  use  the  pi-state  orbitals  to  obtain  a 
ground  state  wavefunction.  In  practical  cases,  two  different  values  of  the 
dipole  transition  moment  are  obtained.  As  the  MC  space  is  enlarged,  increasing 
the  flexibility  of  the  wavefunction,  the  discrepancy  should  decrease  since, 
in  the  limit  that  the  MC  space  becomes  the  full  Cl,  all  orbital  choices  give 
the  same  result.  The  transition  dipole  moments  we  report  have  been  evaluated 
using  the  orbitals  optimized  for  the  upper  state.  Because  the  X-to-pi  dipole 
transition  moments  were  very  small,  no  improved  calculations  were  done. 

On  the  other  hand,  the  transition  dipole  moment  between  the  X  state 

and  the  B  £  state  has  been  shown^>^  to  be  on  the  order  of  unity.  The  B 
u 


state  is  known 

experimentally 

to 

be 

bound 

by 

1.01  eV.  However,  MCSCF 

ca Iculat ions 

for 

the  B  state 

(  in 

the 

basis 

set 

described  above)  using  a 

configuration 

space  analogous 

to 

that 

used 

for 

the  pi  state  calculations 

show  no  binding.  The  wavefunctions  of  ref.  7  for  the  B  state  are  slightly 
better,  but  still  rather  poor  in  both  the  valence  region  and  particularly 
in  the  Rydberg  region,  where  their  basis  set  is  inadequate. 

Our  calculations  on  the  B  state  were  modelled  after  the  selection  of 

(19)  configurations  of  ref.  7.  We  exclude  their  last  two  terms  (by  requiring 

that  2o^  remain  doubly  occupied),  but  include  a  total  of  158  configurations 

resulting  from  expanding  the  valence  orbital  space  to  include  a  set  of  2h 

orbitals.  Rather  than  using  the  CAS-MC  (i.e.,  the  full  Cl  within  this  orbital 

space),  configurations  were  allowed  subject  to  (a)  not  more  than  one  electron 

in  the  2n  orbitals  and  (b)  a  total  of  at  least  two  electrons  in  3o  and  3o  . 

g  u 


y-;.i 


■  •  i 
*•--5 
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3  1 

With  this  MC  space,  the  description  of  binding  (relative  to  the  P  +  D 
asymptote)  was  improved  considerably.  Our  value  (relative  to  the  total  energy 
at  a  separation  of  8  aQ)  was  0.74  eV. 

It  is  well  known  that  avoided  curve  crossings  can  cause  difficulties 

in  the  convergence  of  MCSCF  wavefunctions  because  the  solution  at  a  preceding 

point  on  the  potential  curve  is  likely  to  be  a  relatively  poor  starting  point 

when  the  character  of  the  wavefunction  has  changed  significantly.  Moreover, 

we  were  able  to  obtain  multiple  solutions  for  certain  values  of  the 

internuclear  separation  R.  In  particular,  in  preceding  from  large  values 

of  R  to  smaller  ones,  we  obtained  wavefunctions  in  to  R  =  2.1  a  .  Inside 

o 

of  this,  convergence  problems  were  encountered.  After  much  difficulty, 
solutions  were  obtained  for  the  Rydberg- like  portion  of  the  curve  (as  far 
in  as  R  =  1.8  a^) .  Following  the  solutions  from  small  to  large  R-values, 
we  obtained  a  solution  higher  in  energy  than  the  previous  one  at  R  =  2.281, 
and  at  R  =  2.4  the  solution  returned  to  the  one  previously  found.  I.e., 

between  R  =  2.1  and  2.4  we  have  two  solutions.  The  transition  moment  with 
the  X  state  (using  the  B-state  orbitals)  is  large  for  the  "valence-diabatic" 
wavefunction  and  smaller  for  the  Rydberg-like  one.  The  ground  state  curve 
is  less  well  described  by  the  latter  orbitals  because  of  its  valence  character. 
If  we  choose  to  use  the  B-state  orbitals  from  the  solution  with  the  lower 
total  energy,  we  obtain  a  discontinuous  curve  for  the  transition  moment  (and 
for  the  ground  state  energy).  We  also  know  that  using  X-state-opt imized 
orbitals  to  evaluate  the  transition  moment  will  not  be  accurate  either;  it 
yields  values  similar  to  those  of  the  valence-diabatic  B-state  orbitals.  Thus, 
treatments  which  do  not  account  for  valence-Rydberg  mixing  cannot  obtain 
correct  values  for  the  transition  moment. 

To  improve  the  accuracy  of  the  moments  and  the  agreement  between  the 
values  obtained  with  the  two  orbital  sets,  we  expanded  the  MC  space  (using 
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:he  same  basis  set)  to  310  configurations.  This  expansion  mainly  involved 
idding  single  excitations  to  an  additional  set  of  pi  orbitals;  we  also  allowed 


:he  3o  3o 
g  u 

l-values  of 


pair  to  be  only  singly  occupied. 
1.8-2. 4  a  .  At  R  =  2.4  we 


We  redid 
were  able 


the  calculations  for 
to  obtain  an  MC310 


o 

javefunction  starting  from  the  MC158  orbital  solution.  However,  this  new 
solution  could  not  be  used  for  smaller  R.  We  redid  the  calculations  for 

l  =  1.8  starting  from  the  MC158  orbitals,  and  we  were  able  to  move  to  larger 
1  and  obtain  an  energy  which  differed  (from  the  point  at  R  =  2.4  aQ)  by  only 
J.0065  eV.  The  energy  difference  between  the  lowest  values  obtained  from 
4C310  and  MC158  at  R  =  2.4  a^  is  only  0.12  eV,  so  the  dipole  transition  moments 
for  larger  R  were  calculated  from  the  smaller  MC. 

The  MC310  configuration  space  used  for  the  B  state  also  facilitates 

3  - 

study  of  the  2  (or  E  state)  potential,  at  least  for  R  >  2.1  a^.  This 

curve  comes  in  from  the  same  asymptotic  limit  as  the  B  state,  but  the  potential 

is  purely  repulsive  in  the  valence  region.  Near  R  =  2.75  aQ  it  undergoes 

an  avoided  crossing  with  the  same  Rydberg  state  that  affects  the  B  state 

near  R  =  2.25  a^.  In  this  case,  a  relative  maximum  occurs  in  the  potential 

near  R  =  2.75  a  .  The  curve  then  follows  the  Rydberg  diabat  until  it  must 

o 

assume  valence  character  near  R  =  2.25  a  .  At  this  point  the  magnitude  of 

o 

the  transition  moment  rises  dramatically  to  values  approaching  unity.  Thus 

as  has  been  previously  observed,^  the  dipole  transition  moments  to  the  B 
and  E  states  undergo  a  crossing  where  their  potentials  have  an  avoided 

crossing.  Unfortunately,  difficulties  in  the  MCSCF  calculations  occur  for 
R  <  2.05,  probably  indicating  an  avoided  crossing  with  a  second  Rydberg 


state . 


The  results  of  the  calculations  are  summarized  graphically  in  Figure 
1.  The  solid  line  marked  by  "X"  and  the  dashed  curve  marked  by  dots  are 


84 


MC48  calculations  for  the  ground  and  pi  states,  respectively.  The  absolute 
values  of  the  X  to  pi  dipole  transition  moments  (M)  are  marked  by  dots  in 
the  upper  portion  of  the  figure.  The  two  sets  of  MC158  results  for  the  B 
state  are  indicated  by  the  open  circles,  and  the  dipole  transition  moment 
is  marked  accordingly.  The  MC310  calculations  for  the  B  and  E  states  are 
marked  by  asterisks  and  squares,  respectively.  The  long'dashed  curve  (marked 
with  squares)  is  the  ground  state  MC310  curve  obtained  using  the 
E-state-optimized  orbitals.  The  unmarked  lines  show  the  available  RKR  data 
for  the  X  and  B  states,  with  the  ground  state  minimum  arbitrarily  plotted 
at  -149.8  au.  A  portion  of  the  data  for  Fig.  1  are  given  in  Table  I. 


s2 


A  basis  set  of  54  atomic  orbitals  was  used  for  S2«  It  consisted  of 

the  6s  4p  contracted  set  of  McLean  and  Chandler®  augmented  by  a  set  of  diffuse 

p  functions^  (exponent  =  0.02)  and  a  set  of  d  orbitals  with  exponent  0.12. 

Although  S2  and  O2  have  analogous  excited  state  spectra,  the  potential 

energy  curves  differ  considerably  in  shape.  The  asymptotic  limits  are  closer 

together  in  S2,  but  the  ground  state  is  less  bound  than  O2 .  The  diabatic 

Rydberg  state  is  also  shallower.  As  a  result,  no  avoided  crossing  occurs 
3  3 

in  the  first  li  state  or  the  B  E  state  outside  of  R  =  3.0  an.  (Calculations 
u  u  u 

without  diffuse  functions  have  been  reported  for  these  two  states. ^ 
Calculations  of  Rydberg  character  in  states  of  other  symmetries  have  also 
appeared.1*)  As  in  02,  the  2  E^  state  has  a  local  minimum. 

The  total  energy  curves  and  dipole  transition  moments  (from  the  ground 


state)  for  these  three  states  are  shown  in  fig.  2. 


Because  the  pi  state 


of  greater  importance  in  S2  (for  dynamical  reasons^)  calculations  were 

e  using  the  MC48  analogous  to  that  for  O2  and  using  a  183  term  MC.  The 

ter  orbital  space  was  similar  to  that  for  the  B  state  (analogous  to  the 

B  state  calculations) ,  except  that  60  and  60  orbitals  were  included 

g  u 

itead  of  a  2n  set  and  single  excitations  from  4o  were  allowed.  The  dashed 

u 

rve  with  data  points  marked  by  dots  is  from  the  smaller  treatment,  and 

:  one  with  pluses  is  from  the  MC183.  The  larger  MC  calculations  are  required 

order  to  obtain  a  very  slight  binding  in  this  state. 

The  B  state  calculations  employed  a  162  term  MC.  This  configuration 

ace  was  the  same  as  the  analogy  to  the  MC158  for  O2,  except  that  5o  and 

g 

.  could  be  unoccupied.  The  B  state  data  is  marked  by  open  circles,  and 

e  long-dashed  curve  is  the  ground  state  resulting  from  the  B-state-optimized 
bitals . 

When  comparable  MC's  are  used  for  the  B  and  pi  states,  their  energy 

rves  nearly  coincide  over  a  significant  range  of  R.  This  is  in  agreement 

th  the  experimental  observation^ ,  13  0f  their  interaction. 

3  - 

The  calculations  for  the  2  E  state  used  162  terms  in  the  valence  region 

u 

d  310  terms  in  the  Rydberg  region.  The  data  are  marked  by  triangles  and 
uares,  respectively.  The  local  minimum  occurs  at  a  shorter  distance 

elative  to  the  ground  state  equilibrium  distance)  than  in  O2  because  the 
state  includes  the  minimum  of  the  Rydberg  diabatic  state  (whereas  the  B 
ate  did  in  O2).  The  location  of  this  potential  and  the  root-flipping 

fficulties  encountered  in  the  MCSCF  calculations  for  smaller  R  indicate 
at  interaction  of  the  Rydberg  diabat  with  the  B  state  near  R  =  3.0  au  is 
obable.  We  note  that  the  magnitude  of  the  dipole  transition  moment  between 
is  stale  and  the  ground  state  is  significant  for  values  of  R  near  their 
nima.  We  expect  that  for  values  of  R  <  3.0  a0  the  dipole  transition  moment 
1  the  2  E  state  will  probably  be  larger  than  that  to  the  B  state. 


B6 


RKR  data  for  the  X  and  B  states  are  Indicated  by  solid  unmarked  lines 
dots.  The  X  state  minimum  has  been  arbitrarily  shifted  to  -795.15  au. 
Selected  data  from  Fig.  2  are  given  in  Table  2. 
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APPENDIX  C 


The  Effect  of  an  Avoided  Crossing  on  the 
O2  Schumann-Runge  Photodissociation  Continuum 


Bruce  C.  Garrett 

Lynn  T.  Redmon,  and  Hichael  J.  Redmon 


Chemical  Dynamics  Corporation 
1550  West  Henderson  Road 
Columbus,  Ohio  43220 


Experimental  observation  of  the  O2  Schumann -Runge  photodissociation 
continuum  dates  back  to  the  early  work  of  Landenburg,  Voorhis,  and  Boyce*. 
Since  then,  many  photometric  techniques  have  been  applied  to  quantitative 
measurement  of  the  absorption  spectrum^  and,  more  recently,  an  electron 
energy- loss  spectrum  has  been  measured  for  this  system. 3  The  absorption 
spectrum  is  asymmetric.  It  peaks  near  the  transition  energy  of  8.61  eV  and 
falls  off  rapidly  for  higher  energies.  The  more  recent  experiments^"^  show 
fine  structure  in  the  spectrum. 

3  -  3  - 

The  first  theoretical  calculations  of  the  O2  X  -  B  photodissociation 
continuum  that  were  based  on  full  potential  energy  curves  were  those  by  Jarmain 
and  Nicholls^  and  Bixon,  Rez,  and  Jortner^.  These  studies  used  realistic 
potentials  for  the  bound  part  of  the  B  state,  but  the  repulsive  parts  of 
these  potentials  were  extrapolations.  Their  computed  continuum  absorption 
spectra  were  much  more  symmetric  than  the  observed  spectrum.  In  an  attempt 
to  reproduce  the  rapid  falloff  at  high  transition  energies,  Allison?  proposed 
that  a  part  of  the  B-atate  potential  was  less  repulsive  than  those  previously 
used  and  that  the  dipole  transition  moment  M^g(R)  should  decrease  rapidly 
as  the  internuclear  distance  R  decreased. 

Bixon,  Rez,  and  Jortner^  attributed  the  fine  structure  in  the  absorption 
spectrum  to  continuum  resonances  resulting  from  the  bound  nature  of  the  B-state 
potential.  Jarmain  and  Nicholls^  also  used  a  bound  potential  for  the  B-state 


and  observed  no  structure  in  their  calculated  Franck-Condon  factors.  Huebner 
et.  al.3  interpreted  the  small  peaks  at  low  energies  as  contributions  from 


transitions  to  a  Rydberg  state. 

Since  the  last  theoretical  calculations  on  this  system,  ab  initio 


potential  energy  curves  for  the  B  state  V  and  the  X-to-B  dipole  transition 

D 

moment  have  been  calculated  as  a  function  of  R.®»^  An  avoided  crossing 
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between  the  B  and  E  states  occurs  near  the  X-state  equilibrium  geometry. 
This  causes  V  (R)  to  have  a  shoulder  in  its  repulsive  wall  and  MVD(R)  to 
change  rapidly  in  this  region.  We  find  that  this  shoulder  in  V  (R)  is 

D 

responsible  for  the  rapid  decrease  and  fine  structure  in  the  absorption 
spectrum  at  high  energies. 

In  this  paper  we  present  continuum  oscillator  strengths  df/dE  as  a 
function  of  transition  energy  AE  calculated  using  the  most  accurate  potential 
curves  and  dipole  transition  moments^  available  and  compare  them  with  the 
"apparent"  df/dE  obtained  from  the  experiments  of  Huebner  et.  al .  3  The 
potential  curves  for  the  X  and  B  states  of  O2  are  shown  in  Fig.  1  and  described 
in  reference  10.  The  bound  parts  of  the  potentials  are  from  RKR  data  (see 
reference  2).  The  repulsive  part  of  the  B-state  potential  is  from  ab  initio 
data^  in  to  R  *  1.8  ac  and  is  an  extrapolation  for  smaller  R.  At  R  *  1.8 
a0  the  B-state  potential  energy  is  11.6  eV  above  the  ground-state  vibrational 
energy  in  the  X-state  potential.  The  continuum  oscillator  strength  goes 
to  zero  for  energies  above  10  eV  \  therefore,  the  potential  is  known  over 
a  sufficiently  wide  range  of  R  for  the  present  calculations. 

The  ab  initio  dipole  transition  moment^  for  the  X-to-B  transition  ^^(R) 
is  shown  as  the  long-dashed  curve  in  Fig.  1  and  described  in  Ref.  10.  As 
a  test  of  the  accuracy  of  the  computed  M^g(R)  we  calculated  absolute  oscillator 
strengths  f^,^,,  for  bound  X-to-B  transitions  for  v"  =  0,  1,  and  2  and  v' 

*  0  to  20.  These  are  compared  with  experimental  resultsH»12  in  Table  I. 
The  accurate  calculations  evaluate  f^,^,,  by  numerical  integration  (over  R) 
of  the  product  of  the  initial  and  final  vibrational  wavefunctions  and  MXg(R)» 
and  the  Franck-Condon  approximation  assumes  M  (R)  is  constant  and  equal 
to  its  value  at  the  equilibrium  geometry  of  the  X  state.  In  general,  the 
calculations  employing  the  Franck-Condon  approximation  (FC)  give  oscillator 


strengths  which  are  smaller  than  the  numerically  accurate  ones.  For  v"  = 
0  and  1  the  errors  are  systematic,  both  calculations  underestimate  experiment 
at  low  v'  and  overestimate  at  high  v'.  Therefore,  the  accurate  results  are 
better  at  low  v'  and  the  FC  ones  are  better  at  high  v'.  The  worst  errors 
are  for  v*  *  17;  for  v"  =  0  the  accurate  and  FC  results  are  too  high  by  43% 
and  23%,  respectively,  and  for  v"  *  1  the  errors  are  29%  and  14%.  For  v'<16 
the  accurate  results  are  within  17.5%  for  v"  =  0  and  8.5%  for  v"  ■  1,  and 
the  errors  are  even  smaller  for  v'<10.  For  v"  =  2  the  errors  are  not  as 
systematic  but  the  accurate  results  agree  within  24%  and  the  FC  results  agree 
within  27%.  For  higher  v',  the  vibrational  wavefunction  is  determined  by 
parts  of  the  potential  that  are  less  accurately  known,  the  vibrational 
wavefunction  is  more  diffuse,  and  a  greater  range  of  M^CR)  is  sampled.  The 
overall  agreement  with  experiment  is  acceptable,  especially  at  low  v',  so 
MXB(R)  should  be  sufficiently  accurate  for  calculating  nearly  quantitative 
values  of  continuum  oscillator  strengths. 

To  elucidate  the  effect  of  the  avoided  crossing  in  the  B-state  potential 
on  df/dE,  calculations  were  also  performed  using  a  B-state  potential  curve 
in  which  the  avoided  crossing  was  removed.  This  curve  was  obtained  by 
extrapolating  the  RKR  data  using  the  functional  form  V  (R)  ■  oR  *  exp(-BR), 

D 

where  the  parameters  are  determined  by  matching  the  points  at  the  two  smallest 
R  values  available.  This  model  potential  curve  is  shown  as  a  dashed  line 
in  Fig.  1. 

The  computed  df/dE  are  shown  in  Fig.  2;  the  solid  and  short-dashed  curves 
are  those  obtained  using  the  accurate  and  model  potentials,  respectively. 
These  calculations  are  for  transitions  (of  energy  AE)  from  the  ground 
vibrational  level  (v"  -  0)  of  the  X-state  to  the  B  state;  rotational  motion 
is  neglected.  The  results  using  the  accurate  B-state  potential  show  a  much 
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more  rapid  falloff  at  higher  energies  and  also  exhibit  an  interesting  second 
peak  past  9.0  eV.  The  results  using  the  model  B-state  potential  lack  both 
these  features  and  exhibit  only  the  usual  nearly  symmetric  single  maximum. 

The  lack  or  presence  of  these  features  can  be  better  understood  by 
considering  the  overlap  of  the  ground  vibrational  and  continuum  wavefunctions, 
Xq(R)  and  xe(R)»  respectively.  For  the  model  B-state  potential,  the  major 
peak  in  Xe(R)  shifts  gradually  to  the  left  as  the  transition  energy  is 
increased.  At  energies  near  7.1  eV,  only  the  exponential  tail  of  xE(R) 
overlaps  with  X0(R)  80  df/dE  is  small.  As  the  energy  is  increased  the  overlap 
increases  as  the  major  peak  becomes  centered  over  Xo^R^*  The  overlap  then 
decreases  because  of  negative  contributions  from  the  next  peak  in  xE(R)* 
Near  9.4  eV  df/dE  goes  to  zero  when  the  contributions  perfectly  cancel;  above 
9.4  eV  df/dE  remains  small  because  of  this  cancellation. 

The  maximum  in  df/dE  for  the  accurate  potential  occurs  at  a  higher  energy 
(than  for  the  model  potential)  because  the  accurate  potential  is  more  repulsive 
for  energies  below  8.9  eV.  The  rapid  drop  in  df/dE  and  its  second  maximum 
are  better  understood  by  examining  Fig.  3  which  shows  X£(R)  for  the  transition 
energies  8.9,  9.04  and  9.14  eV,  corresponding  to  the  extrema  in  df/dE.  As 
the  energy  is  increased  from  the  B-state  energy  at  the  avoided  crossing  (8.9 
eV),  Xc.(R)  shifts  rapidly  to  the  left  and  its  major  peak  moves  from  a  position 
of  favorable  overlap  to  one  of  less  favorable  overlap  with  the  ground-state 
vibrational  wavefunction  xq(R)»  accounting  for  the  first  maximum  and  zero 
of  df/dE.  The  second  peak  in  df/dE  is  caused  by  the  overlap  of  the  un-usually 
large  second  peak  in  X-(R)  with  x  (R). 

In  the  previously  described  calculations  of  df/dE,  we  numerically 
integrated  M  (R),  x  ( R)  and  xP(R)  over  R.  In  Fig.  2  we  also  show  the  results 
of  calculations  employing  the  Franck-Condon  approximation  as  the  long-dashed 


curves.  Although  there  are  quantitative  changes  in  df/dE  at  low  energies, 
the  qualitative  features  discussed  above  are  still  present  in  these 
calculations.  This  indicates  that  these  features  are  the  result  of  the  shape 
of  Vg(R)  and  do  not  depend  strongly  upon  the  shape  of  M^dO. 

To  compare  with  the  experimental  results  for  df/dE,  we  average  our 
computed  df/dE  over  a  300  K  and  1000  K  distribution  of  initial  vibrational 
and  rotational  states  and  sum  over  final  orbital  angular  momentum, 


<df/dE>  T  *  £  e‘Gviji/kBT  £  dfVijl*f/dE  (1) 

vUi  *f 

where  the  rovibrational  partition  function  is  given  by 

V  2  .'eV‘1‘/kBT  (2) 

viJi 

eViji  is  the  rovibrational  energy  level  of  the  X  state  for  vibrational  state 

and  rotational  state  J^,  kg  is  Boltzmann's  constant,  and  dfv^j^^^/dE  la 

the  continuum  oscillator  strength  per  unit  energy  for  transition  from  state 

vjJj  of  the  X  state  to  the  B  state  with  orbital  angular  momentum  fcf.  (Note 

that  only  odd  j  states  are  allowed  for  the  X  state  and  Aj**±l  transitions 

are  optically  allowed.)  The  computed  <df/dE>j  for  the  accurate  potential 

3 

is  compared  with  experiment  in  Fig.  4.  Although  the  calculated  peak  is 
a  little  too  high  and  too  narrow,  the  falloff  at  high  energies  is  in 

qualitative  agreement  with  experiment. 

In  conclusion,  the  rapid  decrease  of  df/dE  and  the  presence  of  structure 
at  higher  energies  are  due  to  an  avoided  crossing  in  the  B-state  potential 

curve.  A  rapid  variation  in  M  (R)  is  not  necessary,  as  evidenced  by  the 

AO 

validity  of  the  Franck-Condon  approximation  for  this  transition. 
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i  iVhilfcii.1 


experiment:  for  v"=0  and  1,  M.  Ackerman,  F.  Biaume,  and  G.  Kockarts,  Planet.  Space  Sci.  _18>  1&39  (1970);  for 

R.  D.  Hudson  and  V.  L.  Carter,  J.  Opt.  Soc.  Amer.  _58,  1621  (1968). 


Lpole  transition  moments  vanish  asymptotically.  is  extended  to 

arge  R  values  by  fitting  the  functional  form 


MoiaftR)  =  C  exp(-DR) 


(3) 


o  reproduce  ^^^CR)  at  the  two  largest  R  values  in  the  interpolated  region, 
ubroutines  for  generating  the  transition  dipole  moments  and  potentials 
sed  in  this  work  are  available  from  the  authors. 


02 


Transitions  were  considered  from  the  ground  electronic  state  of  02  (X-^E  ) 

3  - 

:o  three  electronic  states:  the  two  lowest  states  of  symmetry  (labeled 

3 

1  and  E)  and  the  lowest  state  of  symmetry.  The  potential  curves  for 

:hese  states  are  illustrated  in  Fig.  1. 

For  the  X  state,  the  RKR  data  compiled  by  Krupenie^  was  used  in  the  region 

O  O 

From  0.9761  A  to  1.7915  A  and  the  results  of  Saxon  and  Liu2^  were  used  from 

O  O 

1.852  A  to  5.292  A.  The  Saxon-Liu  results  were  shifted  to  have  the 

O 

sxperimental  asymptotic  value  at  10.584  A. 

The  RKR  data  for  the  B  state^  were  used  for  values  of  the  internuclear 

°  °  ?■; 

separation  between  1.33515  A  and  2.57557  A.  The  results  of  Saxon  and  Liu^-" 

O 

*ere  used  from  2.646  to  5.292  A.  The  ab  initio  points  were  shifted  to  have 

O  O 

:he  experimental  asympotic  value  at  10.584  A.  From  0.953  to  1.270  A,  the 

esults  of  Redmon2^  were  used.  The  points  were  adjusted  to  match  the  inter- 

O 

>olated  RKR  data  at  1.376  A. 

No  RKR  data  are  available  for  the  ^Jlu  and  E  states  so  the  fits  were  to 

O 

ab  initio  data.  For  the  ^flu  state,  the  data  points  at  0.953  A  and  for 
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II.  INPUT  DATA 


A.  Potential  energy  curves  and  dipole  transition  matrix  elements 

All  potential  energy  curves  and  transition  dipole  moments  are  fitted 
as  described  in  paper  Is  the  potential  is  fitted  to  the  available 
experimental  and/or  ab  initio  data  points  by  a  cubic  spline  function  and 
represented  at  smaller  R  values  by  the  functional  form 

V(R)  =  A  R*1  exp(-BR)  (1) 


where  the  parameters  are  determined  by  requiring  that  Eq.  (1)  reproduce 
the  two  innermost  data  points.  For  R  values  larger  than  the  spline-fit 
region,  the  potential  has  the  functional  form 


V(R)  =  Vq  -  C6  R'6  -  C8  R-8 


(2) 


where  Vq  is  the  experimental  dissociation  energy  and  the  parameters  are 
determined  by  requiring  Eq.  (2)  to  reproduce  the  outermost  potential  data 
points.  The  cubic  spline  function  is  fitted  so  that  the  first  derivatives 
at  the  left  and  right  ends  of  the  spline-fit  range  match  the  derivatives 
of  the  functional  forms  of  Eqs.  (1)  and  (2),  respectively. 


Over  the 

range  of 

internuclear 

distances 

for 

which  dipole  transition 

moments  are  available, 

fitted 

by 

fourth-order  Lagrange 

.  •» 

interpolation. 

For  smaller  values  of  R,  the 

dipole  transition  moment  is 

assumed  to  be 

constant 

and  equal  to 

the  value 

at 

the  smallest  R  value  for 

which  data  is  available.  For  all  of  the  transitions  considered  here,  the 
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In  this  paper,  electron-impact  cross  sections  for  optically-allowed 
transitions  of  ground  state  O2  and  S2  to  bound  and  dissociative  states  are 
presented.  The  cross  sections  are  for  conditions  relevant  to  laser  plasmas: 
initial  relative  translational  energies  are  taken  from  threshold  to  25  eV, 
initial  vibrational  states  are  selected  from  the  lowest  few  states,  and 
rotational  distributions  are  characterized  by  temperatures  from  0  to  1000 
K. 

The  impact  parameter  (IP)  method  for  diatomic  molecules  as  presented 
by  HaziH’12  i8  designed  to  treat  optically  allowed  transitions.  In  the 
preceding  paper^  (paper  I)  the  present  authors  extended  the  IP  method  to 
treat  excitation  of  vibrationally  and  rotationally  excited  initial  electronic 
states  to  either  dissociative  or  bound  energy  levels  of  the  final  electronic 
states.  The  IP  method  has  been  shown  to  give  reliable  estimates  of  integral 
cross  sections  for  high-energy  collisions  which  are  comparable  to  those 
obtained  from  the  Born  approximation  and  other  "plane-wave"  methods®* 1^”23> 
and  it  agrees  well  with  theoretical  methods  which  are  not  based  upon 
plane-wave  approximations^*  13.  Although  the  assumptions  of  the  IP  model 
may  often  not  be  valid  at  low  energies,  it  can  still  provide  qualitative 

information  about  the  effects  of  vibrational  and  rotational  energy  on  the 
cross  sections  in  this  region. 

The  impact  parameter  method  is  described  in  detail  in  paper  I.  Section 
II  describes  the  data  input  into  the  IP  calculations.  Both  RKR 
(Rydberg-Klein-Rees^)  and  ab  initio  electronic  structure  data  are  combined 
and  fitted  to  obtain  the  best  possible  potential  curves.  The  dipole 
transition  matrix  elements  are  obtained  from  ab  initio  calculations.  The 

Born  calculations  used  to  calibrate  the  IP  calculations  are  also  presented 

in  Sec.  II.  In  Sec.  Ill  the  results  of  the  dynamical  calculations  are 

presented.  Section  IV  discusses  the  results  and  Sec.  V  summarizes  the  study. 


I .  INTRODUCTION 


The  electronic  excitation  and  dissociation  of  O2  in  collisions  with 
low-energy  electrons  is  an  important  process  in  the  photochemistry  of  the 
upper  atmosphere  and  in  some  gas-discharge  lasers.  For  example,  in  gas 
lasers  dissociative  excitation  can  significantly  affect  the  efficiency  of 
the  laser.  Although  the  spectroscopy  of  the  oxygen  molecule  has  been  studied 
extensively,^  there  have  been  only  a  few  experimental  measurements  of  cross 
sections  for  electronic  excitation  by  electron  impact.  Lawrence^  and  Mumma 
and  Zipf^  have  measured  electron-impact  cross  sections  for  dissociation 
through  excited  electronic  states  which  are  higher  than  15  eV  above  the 
ground  state.  Linder  and  Schmidt**  have  measured  integral  cross  sections 

for  spin- forbidden  transitions  to  bound  vibrational  states  and  Trajmar  et. 
al.^»6  have  measured  differential  and  integral  cross  sections  for 

spin- forbidden  transitions  to  bound  vibrational  states.  More  recently, 
Wakiya^  has  measured  differential  and  integral  cross  sections  for  the  X-to-B 
and  some  spin- forbidden  transitions  and  there  has  been  one  calculation  by 
Chung  and  Lin®  of  the  cross  section  for  dissociation  through  the  B  state 
of  O2.  In  the  present  paper  we  are  interested  in  optically  allowed 
transitions  to  bound  vibrational  states  and  to  dissociative  states  of 

low-lying  electronic  states  which  are  accessible  in  collisions  with  low-energy 
electrons  that  are  present  in  gas-discharge  lasers. 

The  sulfur  molecule  is  a  possible  candidate  for  a  gas  laser^  and  has 

been  extensively  studied  spectroscopically,!®  but  little  is  known  about 
its  cross  sections  for  electron-impact  excitation  and  dissociation.  Because 
the  O2  and  S2  molecules  are  isoelectronic ,  it  is  interesting  to  compare 
the  analogous  electron- impact  cross  sections  of  these  systems. 
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1 


The  impact  parameter  method  is  used  to  calculate  integral  cross  sections 

for  electronic  excitation  and  dissociation  of  O2  and  S2  by  electron  impact. 

For  both  molecules,  excitations  to  bound  and  dissociative  states  are 

considered  for  transitions  from  the  ground  electronic  state  (X^E  )  to  the 

g 

two  lowest  states  of  symmetry  (labeled  B  and  E  for  O2  and  B  and  2  for 

S2)  and  the  lowest  state  of  ^IIU  symmetry.  The  dependence  of  the  cross 

sections  upon  initial  vibrational  and  rotational  state  is  studied  for  low 
collision  energies  (threshold  to  25  eV).  For  some  transitions  a  change  in 

initial  vibrational  state  can  have  a  significant  effect  upon  the  cross 

sections,  but,  in  general,  the  effect  of  changing  the  initial  rotational 


state  is  small. 
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Electronic  excitation  and  dissociation  of  0^  and  Sg  by  electron  impact 
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AE  (eV) 


Potential  energy  curves  for  the  X  and  B  states  of  02  and  the  transition 
dipole  moment  for  the  X-to-B  transition.  The  two  solid  curves  are  the 
accurate  fits  to  the  X-  and  B-atate  potentials  from  Ref.  10.  The  short 
dashed  curve  is  a  model  potential  for  the  B  state  which  has  no  avoided 
crossing.  The  long-dashed  curve  is  the  fit  to  the  ab  initio  transition 
dipole  moment  of  Ref.  9.  (An  atomic  unit  of  dipole  moment  is  equal  to 
2.54177  Debye.) 

Continuum  oscillator  strength  per  unit  energy  as  a  function  of  transition 
energy  for  the  X-to-B  transition  in  02.  The  solid  and  long-dashed  curves 
are  obtained  using  the  accurate  B-state  potential  and  the  short-dashed 
curve  uses  the  model  B-state  potential.  The  solid  and  long-dashed  curves 
are  evaluated  by  accurate  numerical  integration  over  R  of  the  initial 
and  final  wavefunctions  and  transition  dipole  moment,  whereas  the 
long-dashed  curves  employs  the  Franck-Condon  approximation  using  the 
transition  dipole  moment  at  the  X-state  equilibrium  geometry. 

Continuum  wavefunctions  for  the  B  state  of  02  at  three  transition  energies: 
a)  8.9  eV,  b)  9.04  eV,  and  c)  9.14  eV.  Also  shown  in  each  part  is 
the  B-state  potential  curve  and  the  ground-state  vibrational  wavefunction 
for  the  X  state.  The  amplitudes  of  the  wavefunctions  are  arbitrary  and 
the  amplitude  of  the  continuum  wavefunctions  have  been  sealed  by  a  factor 
of  0.08  relative  to  the  vibrational  wavefunction. 

Continuum  oscillator  strength  per  unit  energy  as  a  function  of  transitin 
energy  for  the  X-to-B  transition  in  02.  The  solid  curve  is  the 
experimental  data  of  Huebner  et.  al.  (Ref.  3)  and  the  two  dashed  curves 
are  computed  for  rotational  temperature  of  300  and  1000  K. 
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1.376-2.646  A  were  taken  from  Ref.  26  and  the  points  in  the  range  0.995-1.312 

O 

A  were  taken  from  Buenker  and  Peyerimhof f .27«2®  The  ab  initio  points  were 
shifted  to  match  the  experimental  asymptote  and  to  fit  smoothly  from  1.312 

O 

to  1.376  A. 

For  the  E  state,  the  data  points  of  Ref.  26  were  used  in  the  range  from 

O  O 

0.953  A  to  3.175  A  and  the  results  of  Saxon  and  Liu2^  were  used  from  3.440 

O 

to  5.292  A.  The  ab  initio  points  were  shifted  to  match  the  experimental 

O  O 

asymptote  at  10.584  A  and  to  match  each  other  at  3.175  A. 

All  three  excited  state  curves  exhibit  either  a  shoulder  or  a  local  well 

O 

in  the  potential  energy  curves  between  1.111  and  1.164  A.  These  features 
are  the  result  of  avoided  crossings,  where  the  electronic  wavefunction  changes 
character.  For  example,  the  B  state  wavefunction  changes  from  a  state  of 
Rydberg  character  to  the  left  of  the  shoulder  to  one  of  valence  character 
to  the  right  of  the  shoulder. 

The  transition  dipole  moments  used  in  the  cross  section  calculations 

are  plotted  in  Fig.  2.  These  are  all  fits  to  the  ab  initio  results  of  Ref. 

3 

26,  which  are  similar  to  those  of  Ref.  28.  For  the  X  to  B  and  X  to  IT 

u 

O 

transitions,  the  ab  initio  data  was  available  from  0.953  to  2.646  A;  for  the 

O 

X  to  E  transition  the  ab  initio  data  was  available  from  1.111  to  2.249  A.  The 
dipole  transition  moment  changes  rapidly  with  increasing  internuclear  distance 
in  the  region  of  the  avoided  crossing.  This  type  of  behavior  is  not  unexpect¬ 
ed  since  the  electronic  wavefunction  changes  character  rapidly  in  this  region. 


2.  S, 


The  S2  molecule  is  isoelectronic  to  O2,  so  the  same  transitions  were 

3  - 

considered:  from  the  ground  electronic  state  XI  to  the  two  lowest  states 

8 

3  -  3 

of  1^  symmetry  (B  and  2)  and  the  lowest  state  of  ITu  symmetry  (B").  The 
potential  curves  for  these  states  are  illustrated  in  Fig.  3. 
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For  the  X  state,  the  RKR  data  of  Brabson  and  Volkmar^^  was  used  in  the 

O 

range  1.5870  to  2.5456  A.  The  RKR  data^  for  the  B  state  was  U6ed  in  the 

O 

range  1.8093  to  3.0554  A,  and  the  ab  initio  results  of  Ref.  26  were  used 

o 

in  the  range  1.640  to  1.799  A.  The  ab  initio  points  were  shifted  to  match 

0  3  -  3 

an  extrapolation  of  the  RKR  data  at  1.799  A.  The  2  E  and  B"  n  states 

u  u 


used  the  results  of  Ref.  26  in  the  range  1.588  to  3.969  A. 


The  results 


for  the  23E~  were  shifted  to  reproduce  the  experimental  dissociation  energy 
asympotically  and  the  B"  state  results  were  shifted  to  reproduce  the 

O 

experimental  dissociation  energy  at  3.969  A. 

No  evidence  of  avoided  crossings  in  the  B  and  B"  states  was  found  in 

3  - 

the  region  of  the  potential  below  10  eV,  but  the  2  E^  state  exhibits  a  local 
well  above  the  dissociation  limit  much  like  the  E  state  of  02*  The 
equilibrium  geometries  of  the  B  and  X  states  of  S2  are  more  nearly  equal 
than  those  of  O2,  so  the  energy  of  the  S2  B-state  at  the  X-state  equilibrium 
geometry  lies  below  the  B  state  dissociation  limit. 

The  transition  dipole  moments  used  in  the  cross  section  calculations 
are  plotted  in  Fig.  4.  All  were  fits  to  the  ab  initio  results  of  Ref.  26. 
For  the  X  to  B  transitio..,  ab  initio  data  was  available  from  1.693  to  3.969 

O  O 

A;  for  the  X  to  B"  transition,  the  data  extended  from  1.588  to  3.969  A;  and, 

3  - 

for  the  X  to  2  E  ,  transition  the  data  extended  from  1.588  to  2.910  A. 
u 


B.  Minimum  impact  parameter 


In  the  impact  parameter  method,  a  non-zero  bQ  is  used  as  the  lower  limit 

in  the  integration  of  the  transition  probability  over  initial  impact  parameter 

11  13 

to  prevent  divergence  of  the  cross  section.  ’  The  minimum  impact  parameter 
is  determined  by  requiring  that  the  results  of  the  impact  parameter  method 
agree  with  those  of  the  Born  approximation  (BA)  at  high  energies  (800 


eV  in  the  present  application). 


For  comparison  of  the  BA  and  IP  results 


114 


at  very  high  energies  the  following  approximations  are  made:  rotational 

and  vibrational  states  are  treated  as  being  degenerate  and  the  Franck-Condon 
approximation  is  assumed  to  be  valid.  Within  these  approximations,  the 
cross  sections  calculated  are  for  purely  electronic  transitions  summed  over 
all  final  vibrational  and  dissociative  states,  and  the  IP  method  reduces 
to  the  method  as  originally  formulated  by  Hazi^.  At  sufficiently  high 
energy  the  IP  cross  sections  become 


oai“f(E)  - 


2  it  f 


aiaf 


E|  AEotiCtf  1 


In 


b0  l^otjafl 
(2E)iS 


(4) 


where  AEaia  is  the  transition  energy  between  electronic  states,  E  is  the 
electron  translational  energy,  and  is  the  electronic  oscillator 

strength,  which  is  related  to  the  dipole  transition  moment  by 

^laf  (Re>l  (5) 

.  f  gl 

The  sum  is  over  projections  of  the  initial  and  final  electronic  angular  momentum 
along  the  body  fixed  axis,  m  and  e  are  the  mass  and  charge  of  the  electron, 
h  is  Planck's  constant,  and  g£  is  the  degeneracy  of  the  initial  electronic 
state.  In  the  Franck-Condon  approximation  is  assumed  to  be  a  slowly 

varying  function  of  R  and  is  evaluated  at  a  fixed  R  value,  the  equilibrium 
geometry  Re  of  the  ground  electronic  state  in  this  case.  For  a  given 
transition  at  a  fixed  translational  energy  E,  fajOf  an<*  ^a^af  are  known 
and  the  rhs  of  Eq.  (4)  is  equated  to  the  calculated  BA  cross  section  and 
solved  to  obtain  bg. 

The  input  required  for  the  Born-approximation  calculations  consists  of 
the  transition  energy  AEaia^,  the  atomic  basis  functions,  the  transformation 


2m 


ai“f 


3h2e2 


AE, 


A,A 


matrix  for  constructing  the  molecular  orbitals,  and  the  transition  density 


matrix.  The  details  of  the  Born  calculations  are  as  described  by  McCurdy 


and  McKoy  .  21 

The 

transition  energy 

is  obtained  from 

the 

f  its 

to 

the 

potential  curves 

and 

is  the  vertical 

transition  energy 

at 

Re  * 

and 

the 

electronic  structure  information  is  obtained  from  MCSCF  calculations. 26 


The  minimum  impact  parameters  obtained  from  these  calculations  are  presented 
in  Table  I  along  with  the  transition  energies  AEai0(^.  All  bg's  were  evaluated 

O 

at  E  =  800  eV.  For  O2,  all  were  evaluated  at  Re  =  1.207  A  and  for  S2»  values 

O 

of  bo  were  evaluated  at  Re  =  1.889  A. 


III.  RESULTS 


First,  the  vibrational  dependence  of  the  electron  impact  cross  sections 
was  examined.  In  these  calculations  the  rotational  levels  were  assumed 
to  be  degenerate  and  the  cross  sections  for  excitation  to  bound  and 
dissociative  states  were  calculated  from  Eqs.  (31)  and  (40)  of  paper  I, 
respectively.  These  equations  represent  the  IPV  (bound  bound)  and  IPVD 
(bound  -*•  dissociative)  extensions  to  the  original  IP  method. ^  The  effect 
of  rotational  excitation  is  discussed  in  section  III.C. 


A.  02 


Excitation  to  the  BJE  state  of  0o  is  predominately  dissociative.  Figures 

g 


5  and  6  present  the  cross  sections  for  excitation  to  bound  states  and  the 
dissociative  continuum,  respectively.  Dissociation  is  favored  for  this 
state  because  its  potential  in  the  region  of  the  ground  state  equilibrium 
geometry  is  repulsive  and  the  overlap  of  the  resultant  continuum  wavefunction 


with  the  ground-state  vibrational  wavefunction  is  better  than  its  overlap 

with  the  bound  vibrational  energy  levels  of  the  excited  electronic  state. 

As  the  initial  vibrational  energy  is  increased,  the  vibrational  wavefunction 

becomes  more  diffuse  and  its  overlap  with  the  bound-state  wavefunctions 

of  the  excited  electronic  state  improves,  causing  a  large  enhancement  in 

this  cross  section.  Near  threshold,  the  dissociative  cross  sections  also 

increase  with  increasing  initial  vibrational  energy,  mainly  because  the 

energetic  threshold  is  decreased  as  v  ^  is  increased.  However,  they  gradually 

decrease  with  increasing  initial  vibrational  energy  for  higher  energies. 

The  IP  results  for  v^*0  agree  well  with  those  computed  by  Chung  and  Lin® 

using  the  Born-Ochkur  approximation  and  with  the  experiment  results  of  Wakiya.7 

3 

The  cross  sections  for  dissociation  through  the  lowest  Jl^  state  are 
shown  in  Fig.  7.  These  cross  sections  are  typically  twenty  to  thirty  times 
smaller  than  those  for  dissociation  through  the  B  state,  largely  because 
of  the  smaller  transition  dipole  matrix  elements  (see  Fig.  2.) 

The  interpretation  of  the  cross  sections  for  excitation  to  the  second 

3  + 

excited  state  (the  E  state)  is  complicated  by  the  local  well  in  the 

potential  energy  curve  above  the  dissociation  limit  (see  Fig.  1).  The 
potential  barrier  is  1.5  eV  above  the  local  minimum,  and  the  well  supports 
six  quasibound  states,  which  will  appear  as  resonances  in  the  elastic  cross 


section  for 

collision  of 

two 

0  atoms. 

We 

have 

performed 

stabilization^® 

calculations 

to  determine 

the 

locations 

of 

the 

quasibound 

states  of  this 

potential . 

Because  this 

local 

well  in 

the 

E  state  lies 

almost  directly 

above  the  equilibrium  geometry  of  the  X  state,  there  is  a  large  overlap 
between  the  vibrational  wavefunctions  of  the  ground  state  and  the  quasibound 
states.  Therefore,  most  of  the  transitions  from  the  X  state  go  to  these 
quasibound  states  of  the  E  electronic  curve,  and  there  will  be  very  little 


direct  dissociation.  However,  dissociation  from  these  quasibound  states 


can  occur  by  tunneling.  The  contribution  of  these  quasibound  states  to 
the  dissociation  cross  section  is  determined  by  the  competition  between 
|  radiative  relaxation  to  the  X  state  and  tunneling  predissociation. 

The  radiative  lifetimes  of  these  states  have  been  obtained  from  the 

Einstein  A  factors  and  are  calculated  to  be  on  the  order  of  10"^®  seconds. 

I  The  predissociation  tunneling  lifetimes  can  be  obtained  from  the  resonance 

31 

linewidths,  which  are,  in  turn,  estimated  from  the  Gamow  formula.  The 

predissociation  tunneling  lifetimes  are  calculated  to  range  from  0.04  seconds 
I  for  the  lowest  quasibound  state  to  10“^  seconds  for  the  third  quasibound 

state.  Thus  the  radiative  deexcitation  will  be  the  predominant  process 
for  the  three  lowest  quasibound  states,  and  they  will  not  contribute  to 
the  dissociation  cross  section. 

The  cross  sections  for  excitation  to  these  quasibound  vibrational  states 

of  the  E  state  are  presented  in  Fig.  8.  The  resonances  that  occur  because 

i  of  the  three  upper  quasibound  states  are  much  broader  (i.e.,  the  linewidths 

are  larger  and  the  predissociation  lifetimes  are  of  the  magnitude  or  shorter 
than  the  radiative  lifetime)  so  they  should  contribute  to  the  dissociation 
5  cross  sections.  Therefore  contributions  for  final  energies  above  the  first 

three  quasibound  states  were  included  in  the  calculations  of  the  dissociation 
cross  sections  shown  in  Fig.  9.  The  experimental  cross  section  of  Wakiya^ 

for  excitation  to  bound  electronic  states  in  the  energy  range  9.7  to  12.1 

eV  from  v^=0  is  also  shown.  The  state  or  states  responsible  for  the 
transition  in  this  energy  range  is  not  definitely  known,  but  it  has  been 

suggested  that  the  E  state  contributes  a  significant  amount  to  this  cross 
section.  The  experimental  cross  section  is  about  2.5  times  lower  than  the 

IP  results.  The  descrepancies  may  be  the  result  of  a  breakdown  in  the  IP 
method . 


The  cross  sections  for  bound-state  excitation  decrease  with  increasing 
initial  vibrational  state  while  the  dissociative  cross  sections  are  greatly 


enhanced.  The  dissociative  cross  sections  from  the  ground  vibrational  level 
are  ten  to  twenty  times  smaller  for  the  E  state  than  those  for  the  B  state, 
but  for  the  vj-2  state  they  are  only  smaller  by  a  factor  two. 


In  S2  the  B-state  potential  curve  in  the  region  of  the  equilibrium  geometry 
for  the  X  state  lies  below  the  dissociation  limit  of  the  B  state;  therefore, 
excitation  to  bound  states  predominates.  The  B-state  excitation  cross 
sections  are  shown  in  Figs.  10  and  11.  The  cross  sections  for  transitions 
to  bound  states  are  two  to  three  orders  of  magnitude  larger  than  the 

dissociation  cross  sections.  Unlike  the  O2  X-to-B  case,  the  vibrational 
dependence  of  the  bound-state  cross  sections  is  very  weak  and  the  dissociation 
cross  sections  are  greatly  enhanced  as  v^  is  increased.  Again,  the  smaller 
cross  sections  are  caused  by  poor  overlap  between  the  initial  vibrational 
wavefunction  and  the  final  wavefunction  for  nuclear  motion  and  larger  cross 
sections  result  from  the  wavefunctions  for  higher  initial  vibrational  quantum 
numbers  being  more  diffuse  and  overlapping  better  with  the  final  wavefunction. 

The  cross  sections  for  dissociation  through  the  B"  state  are  shown  in 

Fig.  12.  These  cross  sections  display  only  a  weak  dependence  upon  initial 
vibrational  state.  Compared  to  the  cross  sections  for  the  O2  X-to-l^n^ 
transition,  the  S2  cross  sections  have  a  much  earlier  threshold,  rise  to 
an  earlier  peak  (only  2  eV  above  threshold),  then  decrease  more  rapidly 
for  higher  energies. 

3  - 

Excitation  to  the  2  E^  state  of  S2  is  very  similar  to  the  excitation 

3  - 

to  the  E  state  of  O2.  The  2  Eu  state  has  a  local  well  with  its  minimum 
almost  directly  above  the  X  state  minimum;  therefore,  most  of  the 
electron-impact  excitation  to  this  state  is  to  the  nine  quasibound  levels 
of  the  local  well.  As  for  the  O2  E  state,  the  contribution  of  these 

quasibound  states  to  the  dissociation  cross  section  is  determined  by  the 
competition  between  the  predissociation  tunneling  and  radiative  deexcitation. 
The  lifetimes  for  predissociation  tunneling  range  from  10^  s  for  vj*0 


to  6x10  ^  s  for  V£®5.  The  radiative  lifetimes  for  these  states  are  all 
approximately  5x10"^  seconds;  therefore,  excitation  to  the  six  lowest 
quasibound  states  is  considered  nondissociative.  The  cross  sections  for 
excitation  to  the  lowest  six  quasibound  states  are  presented  in  Fig.  13 
and  those  for  dissociation  through  the  2^EU  state  are  shown  in  Fig.  14.  The 
dissociative  cross  sections  are  more  than  four  orders  of  magnitude  smaller 
than  the  cross  sections  for  excitation  to  the  quasibound  states,  but  they 
show  relatively  greater  enhancement  with  increasing  v^. 

C.  Rotational  effects. 

The  effect  upon  the  excitation  cross  sections  of  increasing  the  initial 
rotational  state  was  investigated  for  excitations  from  vibrational  levels 
of  the  X  states  of  C>2  and  S2<  The  cross  sections  for  bound- to-bound  and 
bound- to- continuum  transitions  were  computed  from  Eqs.  (29)  and  (38)  of 
paper  I,  respectively.  The  method  treats  the  diatomic  molecule  as  a  symmetric 
top,  with  spin  conserved.  Thus  there  are  no  transitions  allowed  between 
multiplet  levels.  The  effect  of  changing  initial  rotational  states  is 
expected  to  be  much  smaller  than  the  effect  of  increasing  the  vibrational 
state.  Within  the  IP  method,  only  transitions  in  which  the  rotational  quantum 
number  j  changes  by  ±1  are  allowed.  Therefore  increasing  does  not  change 
the  transition  energy  appreciably  and  the  threshold  energy  for  excitation 
changes  only  slightly.  The  main  effect  of  increasing  is  to  increase 

the  effective  potential  by  the  centrifugal  term  j(j+l)/2pR  .  The  increase 
in  the  potential  is  larger  for  smaller  R,  so  the  potential  well  and  the 
vibrational  wavefunction  for  a  given  level  vj  are  shifted  to  larger 
internuclear  distances  and  changed  in  shape  slightly.  However,  changing 

is  a  small  perturbation  on  the  wavefunction  compared  to  changing  v^. 


Ctoss  sections  have  been  calculated  for  initial  vibrational  state  v^=0 
as  a  function  of  translation  energy  E  and  rotational  temperature  Trot  for 
all  of  the  dissociative  transitions  considered.  These  are  shown  in  Tables 
II  and  III.  The  cross  sections  reported  in  Secs.  III. A  and  III.B,  which 
assume  degeneracy  of  the  rotational  states,  agree  with  those  computed  at 
Trot“0  K  using  Eqs.  (29)  and  (38)  of  paper  I.  Further  calculations  were 

performed  to  j  test  the  rotational  dependence  of  v^=l  and  2  and  the  same 

quantitative  trends  were  observed. 

For  O2  the  dissociation  cross  sections  for  the  X-B  and  X-JI  transitions 
show  very  little  dependence  upon  the  rotational  temperature.  However,  for 
the  X-E  transitions  the  cross  sections  increase  by  almost  a  factor  of  2 

as  the  rotational  temperature  is  increased  from  0  K  to  300  K.  This  is  caused 
by  the  centrifugal  potential  decreasing  the  bound  nature  of  the  E-state 
potential  as  jf  is  increased  and  making  the  effective  potential  more 

repulsive,  thereby  creating  better  overlap  between  the  initial  and  continuum 


wavefunctlons .  Since  the  X-to- 

E  transition 

contributes 

very 

little 

to 

the 

total 

vj=0  dissociative  cross 

section  no 

rotational 

enhancement 

will 

be 

seen , 

but  for  v  £=2  the  X-to-E 

transition 

contributes 

30% 

to  40% 

to 

the 

total  dissociation  cross  section  and  a  slight  enhancement  will  be  seen. 

The  rotational  effects  are  very  similar  for  S2,  the  cross  sections  for 
both  the  X-B  and  X-n  dissociative  transitions  show  very  little  dependence 
upon  rotational  temperature  but  the  cross  section  for  the  X-to-E  dissociative 
transition  shows  a  large  enhancement.  For  S2  the  X-to-2  dissociative  cross 
sections  are  much  smaller  than  those  for  the  X-to-n  transitions  for  v^^O-4; 
however,  a  large  rotational  enhancement  will  be  seen  for  the  total  dissociation 
cross  sections  for  these  v^  (ranging  from  a  factor  of  2  to  a  factor  of  8  for 
v^=0  and  v^*^)  as  the  rotational  temperature  is  increased  from  0  K  to  1000  K. 


D.  Discussion 


In  the  IP  method  it  is  assumed  that  the  contributions  from  large  impact 
parameters  dominate  the  cross  section.  The  minimum  impact  parameter  excludes 
small  b  where  the  transition  probabilities  have  unphysical  divergent  behavior. 
The  extremely  small  value  of  bg  for  the  X-to-B"  transition  in  S2»  therefore, 
reduces  the  confidence  in  the  reliability  of  the  IP  method  for  this  transition 
at  low  to  intermediate  energies. 

To  estimate  the  possible  contributions  to  the  cross  section  from  small 
impact  parameters,  an  alternative  treatment  of  the  small  b  region  was  tested. 
Instead  of  assuming  the  contribution  from  this  region  is  negligible,  the 
contribution  is  assumed  to  be  a  constant  for  this  range  of  impact  parameters 
and  the  cross  section  is  approximated  by 

oIP(E)  =  o}P  +  2v 
b0 

where 

%  -  ’;o  P(VE)  <7> 

~  -  IP 

We  adjusted  b^  to  make  a  (E)  agree  with  the  BA  cross  section  at  800  eV. 

A  measure  of  the  importance  of  the  region  between  b=0  and  b^  is  then  given 


by  the  ratio 

of 

-IP 

°b0 

to 

-ip 

0 

The 

rat  io 

,  -IP  _  IP 

of  a,  to  o  gives 

b0 

an  estimate  of 

the  error 

in 

IP 

0 

(E) 

due 

to 

the 

treatment  of  small 

impact  parameter 

contributions  to  the  cross  section.  The  value  obtained  for  b^  was  generally 
less  than  a  factor  of  two  larger  than  b^,  and  the  errors  were  usually  less 
than  50Z. 
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A  source  of  quantitative  errors  in  the  present  description  of  the  X-to-B" 
transition  in  S2  is  the  potential  energy  curve  used  for  the  B"  state.  The 

O 

B"  curve  of  Fig.  3  does  have  a  slight  shoulder  near  2.5  A  that  is  not  an 
artifact  of  the  fitting  procedure,  but  previous  experimental^, 33  and 
theoretical^  evidence  indicates  that  the  B  state  should  be  slightly  bound. 
Previous  (although  not  definitive)  ab  initio  calculations-^  have  obtained 

O 

a  minimum  in  this  curve  near  2.2  A  at  an  energy  only  0.5  eV  below  the 

dissociation  limit.  In  any  event,  it  is  very  likely  that  the  accessible 
portion  of  the  B"  potential  (near  the  equilibrium  geometry  of  the  X  state) 
does  lie  above  the  dissociation  limit  (of  the  B"  state)  so  that 

electron-impact  dissociation  will  still  dominate  over  excitation  to  bound 
states  of  the  well.  The  shape  of  the  potential  will  change,  thereby  changing 
the  quantitative  values  of  the  dissociation  cross  section;  most  notably 
the  threshold  energy  should  be  shifted  to  lower  energies.  However,  the 
qualitative  prediction  that  the  B"  state  will  be  the  major  route  for 

electron- impact  dissociation  will  remain  the  same. 

One  further  complication  that  arises  in  both  the  O2  and  S2  systems  is 
predissociation  by  nonadiabatic  mechanisms.  In  O2  the  B  state  is  known 
to  be  perturbed  by  the  1  n  state, 35,36  andj  similarly,  the  S2  B  state  is 
perturbed  by  the  B"  state.  In  both  cases  the  coupling  arises  from  spin-orbit 
interactions^, 33, 37, 38_  The  jp  method  does  not  account  for  predissociation 
by  nonadiabatic  mechanisms.  For  the  O2  system,  the  cross  sections  for 

excitation  to  bound  states  of  the  B  state  are  small  and  predissociation 

will  not  contribute  significantly  to  the  dissociation  cross  sections.  For 

S2  the  cross  sections  for  excitation  to  the  bound  states  are  larger  than 

the  dissociation  cross  sections,  and  predissociation  could  make  a 

non-negl igible  contribution  to  the  observed  dissociation  cross  section. 

3  - 

Nonadiabatic  predissociation  could  also  be  important  in  the  E  state  of 


3  -  3 

0£  and  the  2  E^  state  of  S2  since  the  1  state  also  crosses  these  E  states. 
However,  much  less  is  known  about  these  states  (compared  to  the  B  states) 
so  it  is  difficult  to  estimate  the  importance  of  nonadiabatic  predissociation 
for  them. 


V.  SUMMARY 

Electron-impact  cross  sections  have  been  calculated  for  excitation  to 

bound  energy  levels  of  excited  electronic  states  and  for  dissociation  through 

excited  electronic  states  of  O2  and  S2  using  the  impact  parameter  method. 

The  cross  sections  were  the  most  strongly  influenced  by  the  magnitude  of 

the  transition  dipole  moment  and  the  overlap  of  the  initial  vibrational 

wavefunction  with  the  final  wavefunction  for  nuclear  motion.  The  overall 

dissociation  cross  sections  were  often  sensitive  to  the  initial  vibrational 

state,  but  they  were  nearly  independent  of  the  temperature  characterizing 

the  thermal  distribution  of  initial  rotation  states.  The  major  route  for 

electron- impact  dissociation  of  the  O2  ground  state  at  intermediate 

translational  energies  (i.e.,  5  to  20  eV  above  threshold)  via  dipole-allowed 

3  - 

transitions  is  through  the  B  E  state.  The  X-to-B  dissociation  cross  sections 

u 

are  only  weakly  dependent  upon  the  initial  vibrational  state.  Although 

3 

the  cross  section  for  the  X-to-  IT  transition  does  show  significant 

u  & 

enhancement  with  increasing  v^,  its  contribution  to  dissociation  is  small; 
therefore,  the  total  dissociation  cross  sections  show  little  dependence 
upon  initial  vibrational  state. 

The  major  dipole-allowed  contribution  to  the  cross  section  for 

e lect ron- impact  dissociation  of  the  ground  vibrational  level  of  the  X 

3 

state  comes  from  transitions  to  the  B"  II 

u 


state . 


The  B"-8tate  cross  sections 


e  only  slightly  enhanced  with  increasing  V£,  but  the  cross  sections  for 

3  - 

ssociation  through  the  B  state  increase  significantly  and  for  v^=2 

e  similar  in  magnitude  to  the  X-to-B"  dissociation  cross  sections, 
erefore,  for  S^,  the  IP  methods  predict  considerable  vibrational  enhancement 
the  overall  dissociation  cross  section. 
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are  larger  at  high  energy  are  for  ■  0  and  1,  respectively.  The 

long-dashed  curve  is  for  v^  ■  2  and  the  other  solid  and  short-dashed 

curves  are  for  v^  =  3  and  4,  respectively. 

11.  Cross  sections  0°*°^  (E)  in  units  of  cm^  versus  electron  translational 

Vi 

energy  E  for  electron-impact  dissociation  from  the  ground  electronic 
state  of  S2  through  the  B  state  from  initial  vibrational  states  v^ 
“  0-4.  The  cross  sections  are  monotonically  increasing  functions  of 
V£  at  all  energies. 

12.  Cross  sections  (E)  multiplied  by  a  factor  of  10^  in  units  of 

cm^  versus  electron  translational  energy  E  for  electron-impact 
dissociation  from  the  ground  electronic  state  of  S2  through  the  B" 
state  from  initial  vibrational  states  v^  =  0-4  The  solid  and 

short-dashed  curves  which  have  the  highest  threshold  energy  and  one 
smaller  at  high  energy  are  for  v^  ■  0  and  1,  respectively.  The 

long-dashed  curve  is  for  v^  ■  2  and  the  other  solid  and  short-dashed 

curves  are  for  vj  ■  3  and  4,  respectively. 

13.  Cross  sections  0**°^  (E)  multiplied  by  a  factor  of  10^  in  units  of 

cm^  versus  electron  translational  energy  E  for  electron-impact  excitation 

from  the  ground  electronic  state  of  to  quasibound  vibrational  states 
3  - 

of  the  2  E  state  from  initial  vibrational  states  V4  =  0-4. 
u  x 

14.  Cross  sections  (E)  in  units  of  cm^  versus  electron  translational 

vi 

energy  E  for  electron-impact  dissociation  from  the  ground  electronic 

3  - 

state  of  through  2  £u  state  from  initial  vibrational  states  Vj  “ 

0-4.  The  cross  sections  are  monotonically  increasing  functions  of 


v^  at  all  energies.  The  solid  curve  which  is  smaller  at  high  energy 
is  for  v^  “  0.  The  short-dashed  curve  which  is  larger  at  high  energy 
is  for  Vi  *  1.  The  long-dashed  curve  is  for  Vi  ■  2  and  the  other  solid 


and  short-dashed  curves  are  for  Vi  “  3  and  4,  respectively. 


vibrational  states  =  0,  1,  and  2,  respectively.  The  dots  are  the 

Born-Ochkur  results  of  Chung  and  Lin  (Ref.  8)  and  correspond  to  v^®0. 
The  p lu8  sign  is  Wakiya's  experimental  cross  section  for  dissociation 


from  the  vj=0  state  (Ref.  7). 

7.  Cross  sections  c^01^  (E)  multiplied  by  a  factor  of  10^  in  units  of 

cm^  versus  electron  translational  energy  E  for  electron-impact 

3 

dissociation  from  the  ground  electronic  state  of  C>2  through  the  1 
state.  The  solid,  short-dashed,  and  long-dashed  curves  are  for  initial 
vibrational  states  v^  =  0,  1,  and  2,  respectively. 

8.  Cross  sections  o”i0lf  (E)  multiplied  by  a  factor  of  10^  in  units  of 

Vi 

cm2  versus  electron  translational  energy  E  for  electron-impact  excitation 
from  the  ground  electronic  state  of  O2  to  quasibound  vibrational  states 
of  the  E  state.  The  solid,  short-dashed,  and  long-dashed  curves  are 

for  initial  vibrational  states  =  0,  1,  and  2,  respectively.  The 

plus  sign  is  Wakiya's  experimental  cross  section  for  excitation  to 

bound  electronic  states  in  the  energy  region  9.7  to  12.1  eV  (Ref.  7). 

It  is  believed  to  be  dominated  by  optically  allowed  transitions, 

especially  to  the  E  state  (see  Ref.  7). 

9.  Cross  sections  0°^°^  (E)  in  units  of  cm2  versus  electron  translational 

vi 

energy  E  for  electron-impact  dissociation  from  the  ground  electronic 

state  of  O2  through  the  E  state.  The  solid,  short-dashed,  and 
long-dashed  curves  are  for  initial  vibrational  states  v^  =  0,  1,  and 

2,  respectively. 

10.  Cross  sections  0*^°^  (E)  multiplied  by  a  factor  of  10^  in  units  of 

cn2  versus  electron  translational  energy  E  for  electron- impact  excitation 
from  the  ground  electronic  state  of  S2  to  bound  vibrational  states 

of  the  B  state  from  initial  vibrational  states  v^  =  0-4.  The  solid 

and  short-dashed  curves  which  have  the  highest  threshold  energy  and 
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FIGURE  CAPTIONS 

1.  Potential  energy  curves  for  selected  electronic  states  of  02-  The 
fits  of  the  curves  are  described  in  section  II. A. 1. 

2.  Dipole  transition  moments  in  atomic  units  as  a  function  of  internuclear 

3  -  3  - 

distance  R  for  transitions  from  the  X  £  state  of  0o  to  the  B  I  state 

g  z  u 

3  -  3 

(solid  line),  to  the  E  £  state  (long-dashed  line),  and  to  the  1  II 

u  u 

state  (short-dashed  line).  The  curves  are  interpolations  of  ab  initio 
data  as  described  in  section  II. A. 1.  (An  atomic  unit  of  dipole  moment 
is  equal  to  2.54177  Debye.) 

3.  Potential  energy  curves  for  selected  electronic  states  of  S2*  The 
fits  of  the  curves  are  described  in  section  II. A. 2.  The  zero  of  energy 
is  the  minimum  of  the  ground  electronic  state. 

4.  Dipole  transition  moments  in  atomic  units  as  a  function  of  internuclear 

3  -  3  - 

distance  R  for  transitions  from  the  X  £  state  of  S  to  the  B  £  state 

g  2  u 

3  -  3 

(solid  line),  to  the  2  state  (long-dashed  line),  and  to  the  B" 

state  (short-dashed  line).  The  curves  are  interpolations  of  ab  initio 

data  as  described  in  section  II. A. 2.  (An  atomic  unit  of  dipole  moment 

is  equal  to  2.54177  Debye.) 

5.  Cross  sections  oa^a^(E)  in  units  of  cm^  versus  electron  translational 

vi 

energy  E  for  electron- impact  excitation  from  the  ground  electronic 
state  of  O2  to  bound  vibrational  states  of  the  B  state.  The  solid, 
short-dashed ,  and  long-dashed  curves  are  for  initial  vibrational  states 
vj  =  0,  1,  and  2,  respectively. 

6.  Cross  sections  o01*01^  (E)  multiplied  by  a  factor  of  10^  in  units  of 

vi 

2 

cm  versus  electron  translational  energy  E  for  electron-impact 
dissociation  from  the  ground  electronic  state  of  O2  through  the  B  state. 
The  solid,  short-dashed,  and  long-dashed  curves  are  for  initial 
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TABLE  D-II.  Cross  sections  for  dissociation  of  O2  by  electron  impact  for  v^=0. 


E(eV) 

Trot  =  9. 

300 

600 

1000 

X-B 

10 

5. 7 ( -17)a 

5. 7(-17) 

5. 7( - 1 7 ) 

5 . 8( -  17 ) 

15 

7.4(-17) 

7 . 5(-17) 

7.5(— 17) 

7.5(-17) 

20 

8. 2(- 17) 

8-3C-17) 

8.3(-17) 

8.3(-17) 

x-n 

12 

6.2(— 19) 

6  -  4  ( - 19 ) 

6.4(-19) 

6.4(-19) 

15 

7.4(-19) 

7.6(-19) 

7. 7 ( — 19 ) 

7. 7(-19) 

20 

7.2C-19) 

7.4(-19) 

7 . 4 ( — 19 ) 

7.5(-19) 

X-E 

12 

3. 9(-19) 

7 . 9 ( — 19 ) 

8.3(-19) 

8.4(-19) 

15 

8.1(-19) 

1.4(-18) 

1.4(-18) 

1.4(-18) 

20 

8.0(-19) 

1.3(-18) 

1.4(-18) 

1.4(-18) 

a  Numbers  in  parentheses  are  powers  of  10. 


TABLE  D- I .  Minimum  impact 

parameters  and  transition 

energies . 

0 

4^oi  to  j  ^ ^ 

Transition 

bQ  (A) 

02 

X  -*■  b3l’ 

U 

0.921 

8.741 

x  ->  e3e’ 

U 

0.138 

10.315 

x  >  i3n 

U 

0.247 

10.403 

S2 

X  b3z" 

U 

0.974 

4.700 

x  -►  2  V 

u 

2.424 

7.480 

x  -*  B,,3n 

2.28  x  10"16 

6.429 

u 
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The  HC£  molecule  plays  an  important  role  in  laser  plasma  systems  such 
as  the  XeCl  laser.  Dissociation  of  HC&  can  have  a  significant  impact  on 
the  efficiency  of  the  laser.  Stevens  and  Krauss^  concluded  that 
photodissociation  by  UV  radiation  should  not  be  an  important  process  for 
this  system;  however,  to  date  no  consideration  has  been  given  to  the  direct 
dissociation  of  HCi,  by  low-energy  electron  impaot.  In  this  note  we  report 
calculations  of  cross  sections  for  electron- impact  dissociation  of  HCi,. 

The  photodissociation  of  HCi,  by  UV  radiation  is  dominated  by  the 
transition  to  the  A*n  state.  The  A  state  of  HCi.  is  the  lowest  energy  state 
that  is  accessible  from  the  ground  electronic  state  by  an  optical  transition. 
This  state  is  unbound  and  the  transition  from  the  ground  electronic  state 
to  it  has  a  large  oscillator  strength.  As  evidence  of  these  facts  continuous 
absorption  for  this  transition  has  been  observed  experimentally. ^ 

Compared  to  photodissociation,  electron- impact  dissociation  is  more 
complicated  because  of  the  possibility  of  exchange  between  the  incident 
electron  and  electrons  in  the  molecule.  Even  so,  electron-impact  dissociation 
of  HCi  will  be  dominated  by  optically  allowed  transitions  for  electron 
collisions  at  sufficiently  high  energies,  for  which  coulombic  repulsive 
interactions  dominate  over  exchange  interactions.  At  extremely  low  threshold 
energies,  exchange  interactions  become  more  important  and  spin- forbidden 
transitions  to  triplet  states  can  become  an  important  route  for  dissociation 
of  HCi..  However,  for  the  low-energy  collisions  characteristic  of  the  gas-laser 


system,  e lectron- impact  dissociation  of  HCSl  will  proceed  mainly  through  the 
A  state. 

The  presence  of  vibrationally  excited  states  of  HCSl,  under  the  conditions 
of  the  gas  laser,  has  been  indicated  by  both  experiment^  and  kinetic  models 
of  the  laser  plasma. ^  Therefore,  it  is  important  to  assess  the  effect  of 
vibrational  excitation  of  the  HC£  molecule  upon  the  dissociation  cross 
sections . 

The  impact  parameter  method  for  diatomic  molecules^ is  designed  to 
treat  optically  allowed  transitions.  When  extended  to  treat  dissociation 
and  to  include  the  effects  of  vibrational  and  rotational  excitation, ? it 
is  well  suited  for  these  studies.  A  detailed  description  of  the  method  and 
our  computational  procedures  are  presented  in  detail  elsewhere. ^ >® 

As  input  to  the  calculations  we  use  RKR  data  for  the  potential  energy 
curve  for  the  X^E+  state.  9  The  potential  energy  curve  for  the  A^Jl  state 
and  dipole  transition  moment  for  the  X  to  A  transition  are  taken  from  the 
ab  initio  calculations  of  Stevens  and  Krauss.*  These  data  are  shown  in  Fig. 
1.  Also  needed  as  input  to  the  calculations  is  the  minimum  impact  parameter 
bg .  We  computed  electronic  wavefunctions  and  the  components  of  the  transition 
density  matrix  for  the  X  and  A  states  at  the  equilibrium  geometry  of  the 
X  state,  to  be  used  in  Born-approximation  calculations.  A  value  of  the  minimum 
impact  parameter  of  bg  =  2.784  ag  was  obtained  by  requiring  the  original 

impact  parameter  method  of  Hazi^  to  match  the  Born-approximation  calculations 
at  a  translational  energy  of  800  eV.  The  structure  calculations  employed 
a  basis  set  of  46  functions  including  those  appropriate  for  ionic  and  Rydberg 
states.  The  configuration  spaces  for  the  MCSCF  calculations  consisted  of 
67  and  57  terms  for  the  X  and  A  states,  respectively.  The  dipole  transition 
moment  was  also  calculated,  using  the  A  state  optimized  orbitals,  and  the 


ilue  obtained  (0.28507  a.u.)  agreed  well  with  the  value  of  Stevene  and  Krauss. 

As  a  test  of  the  accuracy  of  the  potential  energy  curve  for  the  A  state, 

!  compare  calculated  oscillator  strengths  per  unit  energy  (df/dE)  with 

cperiment.  Experimentally,  continuous  absorption  is  observed  to  begin  at 

>000  cm"  1  and  to  reach  a  maximum  at  65500  cm'^.  Our  calculated  df/dE  peaks 

approximately  73000  cm“^.  If  a  threshold  value  for  df/dE  is  defined  as 

hundredth  of  the  maximum  value,  the  threshold  energy  is  approximately  51000 

The  shift  (above  the  experimental  values)  of  both  the  threshold  energy 

id  the  energy  at  which  the  absorption  is  peaked  indicates  that  the  A-state 

>tential  we  are  using  is  probably  too  repulsive.  The  effect  of  this  on 

le  dissociation  cross  section  is  a  shift  to  higher  energy,  but  the  magnitudes 

E  the  cross  sections  should  be  essentially  correct. 

The  cross  sections  for  dissociation  of  HCH  through  the  A  state  are  shown 

l  Fig.  2  for  v^=  0,  1,  and  2.  As  v^  is  increased,  the  threshold  is  moved 

3  lower  energies  and  there  is  an  enhancement  in  the  cross  sections.  Figure 

shows  the  dependence  of  the  cross  section  upon  initial  rotational  state. 

le  cross  sections  in  Fig.  2  correspond  to  transitions  from  the  ground 

Jtational  state.  In  Fig.  3  we  present  cross  sections  averaged  over  Boltzmann 

Lstributions  of  rotational  states  characterized  by  temperatures  of  0,  300, 

)0,  and  1000  K.  These  cross  sections  are  a  monotonically  decreasing  function 

E  T  ;  as  T  is  increased  from  0  to  300  K  the  cross  sections  decrease 
ret  rot 

r  about  a  factor  of  2  but  further  increases  in  T  decrease  the  cross 

rot 

•ctions  by  less  than  5%.  Since  these  cross  sections  are  approximately  an 
rder  of  magnitude  larger  than  the  cross  section  for  photodissociation 
ilculated  by  Stevens  and  Krauss,  e lectron- impac t  dissociation  will  play 
more  significant  role  in  laser  system  containing  HCi. 
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FIGURE  CAPTIONS 


1.  Potential  energy  curves  (solid  lines)  for  the  X  and  B  states  and  (dashed 
lines)  A  and  C  states  of  HgBr. 


2. 

Dipole  matrix  elements  for 

the  X-B 

transition 

(solid  curve), 

and 

X-A 

and  X-C  transitions  (dashed 

curves  ) 

in  HgBr. 

The  curves  are 

fits 

to 

the  data  of  Wadt  (see  text). 

3.  Cross  sections  multiplied  by  a  factor  of  10^  in  units  of  cm^  versus 

electron  translational  energy  for  electron  impact  dissociation  of  HgBr 
through  the  A  state.  The  solid  and  dashed  curves  are  for  initial 
vibrational  states  =  0  and  4,  respectively . 

4.  Cross  sections  multiplied  by  a  factor  of  10^  in  units  of  cm^  versus 

electron  translational  energy  for  electron- impact  dissociation  through 
the  A  state  of  the  =  4  vibrational  state  of  the  X  state  of  HgBr.  The 
solid  and  dashed  curves  are  for  rotational  temperatures  of  0  and  1000 
K,  respectively. 

5.  Cross  sections  multiplied  by  a  factor  of  10*®  in  units  of  cm^  versus 

electron  translational  energy  for  electron  impact  excitation  of  the  B 
state  of  HgBr.  The  solid  and  dashed  curves  are  for  initial  vibrational 
states  =  0  and  10,  respectively. 

6.  Cross  sections  multiplied  by  a  factor  of  10*^  in  units  of  cm^  versus 

electron  translational  energy  for  electron  impact  excitation  of  the  C 
state  of  HgBr.  The  solid  and  dashed  curves  are  for  initial  vibrational 
states  V|  =  0  and  10,  respectively. 
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excitation  of  the  B  state  for  initial  vibrational  levels  *  0  and  10.  There 
is  at  most  a  40%  enhancement  for  vj  =  10  over  =  0. 

In  Fig.  6  bound  excitation  cross  sections  are  shown  for  the  X-C  transition, 
again  for  vj  =  0  and  10.  The  enhancement  near  threshold  is  caused  by  the 
different  rate  of  change  with  energy  of  the  structural  factor  (influenced 
by  Franck-Condon  factors)  and  the  dynamical  factor. ^ 
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The  dipole  matrix 


and  5.3  A,  and  by  extrapolation  outside  of  this  range.® 
elements  are  shown  in  Fig.  2.  These  are  4-point  Lagrange  interpolations  of 
Wadt's  data. 

The  A  state  is  totally  repulsive  and  is  expected  to  be  a  major  contrib¬ 
utor  to  the  overall  HgBr(X)  dissociation  cross  section.  In  Fig.  3  we  present 
cross  sections  for  dissociation  through  the  A  state  from  threshold  to  25  eV 
for  initial  vibrational  states  equal  to  0  and  4.  The  cross  section  rises 
rapidly  to  a  value  near  1.5x10"^®  cm^  at  2.0  eV  and  subsequently  falls  off 
smoothly  with  increasing  energy.  There  is  only  a  small  effect  on  the  magnitude 

of  the  cross  section  caused  by  varying  initial  vibrational  energy,  although 

there  is  some  enhancement  predicted  by  the  impact  parameter  model  near 
threshold. 

In  Fig.  4  the  effect  of  initial  rotational  temperature  is  presented.  The 
rotational  effect  was  examined  for  initial  vibrational  states  from  *  0 

to  4.  The  same  quantitative  trends  were  observed  for  all  and  we  report 
only  those  for  v  ^  "  4.  The  decrease  of  a  factor  of  two  observed  in  going 
from  Trot  =  0  to  1000  K  is  somewhat  of  an  anomaly  compared  to  what  we  have 
seen  in  other  systems^-^  and  results  from  the  fact  that,  for  Trot  ■  °»  only 
j^  =  0  contributes  to  the  cross  section,  while  for  Trot>0,  other  angular  mo¬ 
mentum  states  contribute. ^  The  contributions  from  various  angular  momentum 
components  are  weighted  by  different  symmetry  factors  (for  E  to  n  transitions 
only)  which  results  in  a  low  temperature  enhancement  of  the  cross  section. 

For  the  X-B  and  X-C  transitions  from  initial  vibrational  states  0  to  4 
the  dissociation  cross  sections  are  extremely  small,  being  less  than  10"^® 
cm^  for  the  X-B  transition.  This  is  reflected  in  the  fact  that  the  sum  of 
Franck-Condon  factors  for  transitions  to  bound  vibrational  levels  of  the  B 

and  C  states  is  unity  to  within  the  accuracy  of  the  calculations  (approximately 
six  significant  figures).  In  Fig.  5  we  present  cross  sections  for  bound 


■mi  ■«  "«  ji  r> 


formulation  is  presently  restricted  to  optically-allowed  transitions.  However, 
vibrational  and  rotational  degrees  of  freedom  can  be  treated  explicitly. ? 

The  impact-parameter  method  requires  as  input  potential  energy  curves 
and  dipole  transition  matrix  elements.  In  addition,  a  value  for  the  lower 
bound  of  an  integral  over  the  impact  parameter  (the  minimum  impact  parameter 
bg)  is  required. ^  This  is  obtained  by  requiring  the  impact-parameter  cross 
section  to  equal  the  cross  section  obtained  from  the  Born  approximation^ 
at  high  energies.  For  the  Born  calculation  we  require  the  one-electron  density 
matrix  from  the  electronic  structure  calculation. 

For  HgBr  we  used  the  potential  curves  and  dipole  matrix  elements  from 
the  ab  initio  calculations  of  Wadt.^  These  calculations  employed  relativistic 
pseudopotentials^2  within  the  configuration  interaction  (Cl)  method,  and  includ¬ 
ed  spin-orbit  effects  through  a  semiempir leal  scheme.  The  basis  set  was  of 
double-zeta  plus  polarization  quality.  To  obtain  the  one-electron  density 
matrices  required  by  the  Born  approximation  calculation  (to  fix  bg)  we  performed 
multiconfiguration  self-consistent-field  (MCSCF)  calculations  using  Wadt's 
basis  set  and  pseudopotential  parameters.  The  one-electron  integrals  were 
obtained  with  the  new  integral  package  of  Pitzer.l-*  MCSCF  calculations  were 
performed  for  all  four  states  at  the  experimental  equilibrium  geometry  of 
the  X  state  using  various  configuration  spaces  having  from  one  to  264  terms. 
We  had  difficulty  obtaining  a  single  set  of  orbitals  that  described  the  ground 
and  excited  states  equally  well,  so  we  chose  the  set  that  best  reproduced 
the  transition  dipole  matrix  elements  of  Wadt^  for  these  transitions.  This 
set  was  used  to  determine  the  density  matrices  for  the  Born  calculation.  The 
computed  minimum  impact  parameters  for  the  X-to-A,  X-to-B,  and  X-tu-C  transi¬ 
tions  were  0.64  ag,  2.1  ag,  and  1.9  ag,  respectively. 

The  potential  energy  curves  used  in  the  calculations  are  shown  in  Fig. 
1.  They  are  represented  by  spline  interpolation  of  Wadt's  data  between  2.1 
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Electron  impact  dissociation  of  HgBr  is  investigated  using 
a  modified  impact-parameter  method.  Cross  sections  are 
calculated  for  dissociation  through  the  B(2E+),  A(2JI+) 
and  C(2£+)  states.  The  X-A  transition  is  found  to  dominate 
the  dissociation  process.  There  is  only  a  small  enhancement 
of  the  cross  section  due  to  molecular  vibrational  energy, 
and  no  rotational  enhancement.  Cross  sections  are  also 
presented  for  nondissociative  excitation. 

The  mercury  bromide  (B-->X)  electronic  dissociation  laser  has  recently 
received  considerable  attention  because  of  its  high  efficiency,  scalability, 
and  blue/green  wavelength  (502  nm) .  *  The  B  state  of  HgBr  can  be  obtained 
from  HgBr2  by  photodissociation, 2  by  discharge  techniques,-*  or  by  electron 
beam  pumping.^  A  recent  review  of  the  properties  of  e-beam  controlled  HgBr 
discharges  is  provided  by  Nighan,  along  with  extensive  references.^  In  this 
report  we  are  concerned  with  electron  impact  excitation  and  dissociation  proc¬ 
esses  for  the  X  to  B(2E+),  A(2I1+)  and  C(2I+)  states  of  HgBr  because  of  their 
possible  importance  to  closed-flow,  repet  it ive- pulse  HgBr  laser  kinetics. 

We  have  applied  a  modified  impact-parameter  method**’2  that  is  well  suited 
for  the  calculation  of  post-threshold  electron  impact  cross  sections  in  diatomic 
molecules.  The  method  has  previously  been  applied  to  the  study  of  electronic 
excitation  and  dissociation  processes  in  O2  and  S2,®  and  in  HCl.^  In  this 
method,  the  cross  section  is  expressed  as  the  product  of  an  electronic  structure 
factor  and  a  dynamical  factor.  The  former  is  obtained  from  an  ab  initio  elec¬ 
tronic  structure  calculation  and  the  latter,  representing  the  dynamics  of 
the  scattered  electron,  is  expressable  in  terms  of  simple  functions. ^ 
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Figure  Captions 


1.  Potential  energy  curves  (solid  lines)  for  the  X  and  A  states  of  HC£  and 

the  dipole  transition  moment  (dashed  line,  in  atomic  units)  for  the 

X  to  A  transition. 

2.  Cross  sections  multiplied  by  a  factor  of  10^®  in  units  of  cm^  versus 

electron  translational  energy  for  electron- impact  dissociation  of  HCfc 

through  the  A  state.  The  solid,  short-dashed,  and  long-dashed  curves 
are  for  initial  vibrational  states  v^  =  0,  1,  and  2,  respectively. 

3.  Cross  sections  multiplied  by  a  factor  of  10^®  in  units  of  cm^  versus 

electron  trans lat iona 1  energy  for  electron-impact  dissociation  of  HCH 

through  the  A  state.  The  highest  solid  curve  is  for  a  rotational 
temperature  of  0  K,  the  next  solid  curves,  short-dashed  and  long  dashed 
curves  are  for  300,  600,  and  1000  K  respectively.  All  cross  sections 

are  for  v^=0. 
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APPENDIX  G 


Program  Flowcharts 


Impact  Parameter  Method 


Program  IP 


Description:  Compute  cross  sections  for  electronic  excitation  of  diatomic 
molecules  by  electron  impact  using  the  impact  parameter  method  of  Hazi. 
Vibrations  and  rotations  are  treated  as  being  degenerate  and  the  fixed-nuclei 
approximation  and  Franck-Condon  approximations  are  used.  Cross  sections 
correspond  to  those  summed  over  final  rovtbrational  states. 
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Program  IPV 


Description:  Compute  cross  sections  for  electronic  excitation  of  diatomic 

molecules  to  bound  vibrational  states  by  electron  impact  using  the  impact 
parameter  method  .  Vibrations  are  treated  explicitly  within  the 

fixed-nuclei  approximation  and  rotations  are  treated  as  being  degenerate. 
Cross  sections  are  summed  over  final  vibrational  states,  dissociative  states 
are  excluded. 


Flowchart : 
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Program  IPV  (INPUT  from  data  file) 


Flowchart : 
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INITIALIZE 
potential  and 
constants 


Program  IPV  (Obtain  x 


Flowchart : 


rara  IPVR 


ription:  Compute  cross  sections  for  electronic  excitation 

:ules  to  bound  vibrational  states  by  electron  impact  using 
neter  method  of  Hazi.  Vibrations  and  rotations  are  treated 
in  the  fixed-nuclei  approximation.  Cross  sections  are  summed 
arational  states,  dissociative  states  are  excluded. 
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